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EXECUTIVE  SUMMARY 


Geotechnical  surveys  have  been  recently  completed  in  the  north  boundary 
vicinity  of  Rocky  Mountain  Arsenal  (RMA)  to  define  the  Arsenal's  subsurface 
hydrogeology.  Water  table  measurements  and  permeability  estimates  have  been 
generated  to  refine  flow  direction  and  quantify  groundwater  movement  of  the 
alluvial  aquifer  as  it  flows  across  the  north  boundary.  This  document  pre¬ 
sents  a  compilation  of  data  assessments  based  upon  the  aforementioned  survey 
results  on  expected  fluoride  contaminant  loading  within  the  proposed  expanded 
north  boundary  control  system. 

Groundwater  at  RMA  acts  as  one  continuous  hydrogeologically  connected  unit 
flowing  from  south  to  north.  Locally  in  the  north  boundary  vicinity  the 
alluvial  aquifer  can  best  be  envisioned  as  two  separate  subsurface  water 
units  which  ultimately  make  up  one  alluvial  flow  crossing  the  Arsenal's  north 
boundary.  The  most  potentially  troublesome  groundwater  unit  moves  beneath 
Basin  F  in  a  northeastemly  direction  toward  the  north  boundary/  Contaminants 
leached  from  surface  waste  basins  move  'within  this  subsurface  flow  and  cross 
the  boundary  in  the  immediate  area  of  the  present  pilot  containment/ treatment 
system.  The  other  groundwater  unit  flows  northwesternly  beneath  and  parallel 

10  the  First  Creek  surface  stream.  This  alluvial  pathway  is  relatively  free 
I  contaminants  and  has  a  much  larger  volume  than  does  the  contaminated 
roundwater  flow.  As  both  subsurface  flows  approach  the  north  boundary, 
water  table  contours  straighten  and  become  parallel  to  the  Arsenal  boundary. 

In  this  area  groundwater  flow  is  directly  northward  with  some  contaminants 
crossing  the  boundary  throughout  the  alluvial  aquifer. 

Based  upon  results  of  these  surveys ,  a  decision  has  been  made  by  the  Army 
to  expand  the  1500  foot  north  boundary  pilot  containment/treatment  system. 
Applicable  water  quality  guidelines  (particularly  the  0.2  ppb  State  of 
Colorado  limit  for  DBCP)  dictate  that  the  extension  must  be  across  the  entire 
northward  flowing  alluvial  aquifer.  Thus,  current  plans  call  for  the  present 
pilot  containment  system  to  be  expanded  approximately  3500  feet  eastward  and 
700  feet  westward.  \ 


Interception  of  the  entire  alluvial  aquifer  will  result  in  compositing  both 
the  contaminated  and  relatively  noncontaminated  groundwater  flow  units  described 
above.  For  design  purposes,  assessments  have  been  completed  to  identify  expected 
contaminant  loadings  within  the  composite  stream  to  be  handled  by  the  expanded 
control  system.  In  July  1979  D'Appolonia  Consulting  Engineers  provided  a 
conceptual  design  of  the  extension  of  the  pilot  containment  subsystem  to  the 
east  and  west.  Fluoride  was  among  four  pollutants  examined  to  predict  total 
quantities  of  contaminants  to  be  intercepted.  A  second  assessment,  performed 
by  US  Army  Toxic  and  Hazardous  Materials  Agency  (USATHAMA)  utilized  the  same 
geotechnical  survey  data  base  to  predict  contaminant  loadings  to  the  expanded 
ilot  treatment  subsystem. 
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As  a  result  of  these  assessments,  it  has  been  determined  that  there  is  a 
good  possibility  that  fluoride  levels,  within  the  expanded  treatment  sub-=- 
system  will  be  naturally  below  the  State  of  Colorado  limit.  This  may  negate 
the  need  to  add  a  costly  fluoride  removal  module  to  the  proposed  granular 
carbon  organic  removal  process.  The  results  of  these  assessments  are  summarized 
below: 


Alluvial  Flow  Rate  (gpd) 


Assessment  A; 
D’Appolonia 
Consulting 
Engineers 


Assessment  B; 

US  Army  Toxic 
and  Hazardous 
Materials  Agency 


638393 


882200 


Total  Fluoride  Migrating 

off  N.  Boundary  of  RMA  (kg/day) 


6.6 


7.8 


Expected  Manifolded 

Fluoride  Concentration  (mg/Z) 


2.3 


State  of  Colorado  Standard  (mg/S-) 


2. A 


2.4 


Variations  noted  in  the  estimates  are  due  to  differences  in  the  choice  of 
aquifer  hydrodynamic  parameters.  Permeability  estimates  for  the  most  per¬ 
meable  aquifer  material  range  from  400  to  600  feet  per  day  for  the  D'Appolonia 
and  USATHAMA  assessments,  respectively.  Equal  variation  is  noted  in  saturated 
thickness  estimates.  These  differences  are  within  an  acceptable  range,  however. 
Geohydrologic  definition  is  not  an  exact  science  and  is  commonly  assumed 
adequate  if  an  80  percent  accuracy  is  achieved.  It  is  emphasized  that  both 
assessments  are  considered  valid  within  the  limits  of  the  accuracy  of  the  data. 
Refinement  of  the  flow  and  fluoride  concentration  expectations  will  not  be 
possible  until  actual  expanded  system  operation  is  accomplished. 
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1.0  INTRODUCTION 


The  Rocky  Mountain  Arsenal  (RMA)  is  located  approximately  10  miles 
northeast  of  the  central  business  districfe^of  Denver,  Colorado  and 
immediately  north  of  the  Stapleton  International  Airport  (Figure  1). 

RMA  was  established  in  1942  and  historically  has  either  produced  toxic 
chemicals  and  chemical  filled  munitions,  or  demilitarized  these  same 
items.  In  1946,  a  large  portion  of  the  manufacturing  facilities  was 
leased  to  private  industry  for  the  production  of  herbicides  and  insec¬ 
ticides.  Chemical  wastes  generated  collectively  by  these  operations 
have  been  discharged  into  several  waste  storage  basins  located  on  the 
Arsenal  grounds. 

The  first  reported  indication  of  off-post  contamination  occurred  in  the 

•er  of  1951,  when  some  crop  damage  was  reported  on  an  irrigated 

northwest  of  the  RMA  (Kolmer  &  Anderson,_ 1977) .  In  1954,  several 
farmers  north  of  the  arsenal  complained  of  damage  to  crops  irrigated 
with  water  pumped  from  the  alluvial  aquifer.  Due  to  these  complaints 
and  subsequent  damage  claims,  the  Department  of  the  Army  initiated 
several  studies.  These  studies  resulted  in  the  construction  of  a  new 
disposal  basin  with  a  low  permeability  liner  (Reservoir  "F",  see  Figure 
1).  Since  1957,  .all  chemical  wastes  have  been  pumped  into  this  reservoir. 

In  May  1974,  diisopropyl  methylphosphonate  (DIMP)  and  dicyclopentadiene 
(DCPD)  were  detected  in  waters  discharging  from  a  bog  located  along  the 
north  boundary  of  the  RMA.  DIMP  was  also  detected  in  water  supply  wells 
for  the  city  of  Brighton  in  December  of  1974.  DIMP  is  a  persistent 
compound  .produced  in  small  quantities  during  the  manufacture  of  GB,  a 
chemical  warfare  agent.  DCPD  is  a  chemical  used  in  the  production  of 
insecticides.  The  off-post  detection  of  DIMP  and  DCPD  prompted  the 
Colorado  Department  of  Health  to  issue  three  Cease  Orders  on  April  7, 
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1975,  that  required  an  immediate  stop  to  surface  and  subsurface  discharge 
of  DIMP  and  DCPD,  development  of  a  plan  to  preclude  future  discharge  of 
the  contaminents ,  and  development  of  a  monitoring  program  to  verify  com¬ 
pliance  with  these  orders. 

In  the  summer  of  1976  analysis  of  groundwater  from  the  north  boundary 
also  revealed  the  presence  of  inorganic  fluorides  and  three  organic 
sulfur  compounds  (p— chlorophenyl  methyl  sulfide,  p— chlorophenyl  methyl 
sulfoxide,  and  p— chlorophenyl  methyl  sulfone).  In  1978,  dibromochloro— 
propane  (DBCP  or  Nemagon)  was  discovered  in  the  groundwater  in  the 
vicinity  of  the  north  boundary  of  the  Arsenal.  Although  these  compounds 
were  not  cited  in  the  Cease  and  Desist  Orders,  they  are  included  in  the 
list  of  compounds  requiring  treatment. 

From  1975  to  1977,  several  investigators  were  involved  in  hydrologic 
investigations  and  the  design  of  a  containment  and  treatment 
system  for  a  portion  of  the  northern  boundary  of  the  RMA.  These  studies 
and  reviews  were  conducted  by  Konikow  (1975),  Re3rnolds  (1975),  Miller 
(1977),  Mitchell  (1976),  Kolmer  and  Anderson  (1977a  and  b) ,  Thomas,  et 
al.,  (1977),  and  Robson  (1976).  The  studies  resulted  in  the  installation 
of  the  present  pilot  containment  system  along  a  portion  of  the  northern 
RMA  boundary  (Figure  1). 

The  selected  system  design  for  containment  and  treatment  of  contaminated 
groundwater  relics  on  the  use  of  an  impermeable  barrier  (bentonite 
slurry  wall)  to  impede  the  natural  subsurface  flow  of  ground  water 
across  the  boundary.  Groundwater  flowing  toward  the  barrier  is  removed 
from  the  upgradient  side  of  the  barrier  by  dewatering  wells  and  treated 
for  the  removal  of  organic  contaminants.  The  treated  water  is  then 
injected  into  the  aquifer  on  the  downgradient  side  of  the  slurry  wall  in 
a  line  of  recharge  wells.  A  schematic  representation  of  the  contain¬ 
ment  system  is  provided  in  Figure  2. 
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Since  installation  of  the  pilot  containment /treatment  system,  additional 
geohydrologic  data  has  been  collected  in  the  north  boundary  vicinity. 

This  information  has  been  used  to: 

o  Provide  a  detailed  definition  of  the  geology  and  groundwater 
in  the  alluvium  between  Basin  F  and  the  north  boundary 
(Zebell,  1979). 

o  Provide  a  description  of  the  locations,  movements,  and 
concentrations  of  various  pollutants  in  the  alluvial 
aquifer  (Thompson,  1979) . 

o  Provide  an  evaluation  of  the  north  boundary  pilot  system 

after  1  year  of  operation  (D'Appolonia  Consulting  Engineers, 
1979). 


o  Provide  a  conceptual  design  of  the  extension  of  the  pilot 
containment  system  to  the  east  and  west  (D’Appolonia 
Consulting  Engineers,  1979). 

#  ed  on  the  success  of  the  pilot  system,  a  decision  has  been  made  by  the 
Department  of  the  Army  to  extend  the  present  pilot  containment  treatment 
system.  The  north  boundary  containment  system  expansion  concept  will 
be  designed  to  recover  water  flowing  in  the  alluvial  aquifer  across 
the  north  boundary  of  the  RMA  with  a  minimum  of  changes  in  the  presystem 
head  distribution.  Predictions  indicate  that  an  average  of  450  to 
650  gallons  per  minute  (gpm)  of  groundwater  crosses  the  north  boundary 
under  present  conditions.  The  requirement  to  treat  all  groundwater 
containing  DBCP  concentrations  above  the  detection  limit  of  0.2  ppb 
dictates  that  this  entire  northward  flowing  aquifer  must  be  controlled 
and  intercepted.  A  3500  foot  eastward  and  700  foot  westward  expansion 
to  the  pilot  system  is  currently  planned. 


Particular  attention  has  been  given  to  the  prediction  of  DIMP ,  DBCP, 
and  fluoride  contaminant  loadings  to  the  treatment  system  once  the 
groundwater  is  manifolded  after  dewatering.  This  analysis  is  limited 
the  issue  of  fluoride.  Independent  assessments  have  been  recently 
completed  by  private  contractor  and  Government  personnel.  The  following 
sections  delineate  results  of  these  efforts. 

o  Section  2;  Site  Hydrology. 

o  Section  3;  Assessment  A  performed  by  D'Appolonia  Consulting 
Engineering 

o  Section  4:  Assessment  B  performed  by  US  Army  Toxic  and 
Hazardous  Materials  Agency. 

o  Section  5;  Conclusions 
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2.0  SITE  HYDROGEOLOGY 


The  hydrogeologic  system  of  concern  along  the  north  boundary  of  the  RMA 
consists  of  an  unconsolidated  alluvial  sand  and  gravel  aquifer  that 
overlies  a  much  lower  permeability  shale  and  claystone  bedrock.  Subsur¬ 
face  flow  of  contaminated  groundwater  to'  the  north  takes  place  within 
this  alluvial  aquifer  and  results  in  a  discharge  across  the  north 
boundary  of  the  Arsenal. 


2.1  BEDROCK  COMPOSITION  AND  TOPOGRAPHY 

The  alluvial  aquifer  is  underlain  by  predominantly  shale  and  claystone 
bedrock  of  the  Denver  Formation.  Previous  studies  of  the  groundwater 
contamination  have  assumed  that  the  major  portion  of  flow  takes  place 
within  the  alluvium  due  to  the  extreme  permeability  contrast  between  the 
bedrock  and  unconsolidated  alluvial  units.  A  number  of  deep  borings 
show  that  the  bedrock  is  composed  primarily  of  shale  and  claystone  with 
occasional  silt  lenses.  Some  poorly  cemented  sandstones  have  also  been 


feet. 


including  one  continuous  unit,  about  4  feet  thick,  that  can  be 
through  several  borings  over  a  distance  of  greater  than  2,000 
Other  sandstones  up  to  12  feet  thick  appear  to  be  lenticular. 


A  weathered  zone  is  found  in  the  shales  extending  to  25  to  35  feet 
below  ground  surface  in  the  lower  lying  areas  and  up  to  50  feet  below 
ground  surface  in  the  higher  areas.  This  weathered  zone  generally 
extends  from  5  td  15  feet  below  the  bedrock  alluviiim  contact.  At  the 
eastern  end  of  the  existing  slurry  wall,  the  weathered  zone  extends 
approximately  7  feet  below  the  bottom  of  the  wall. 

The  sandstones  within  the  bedrock  and  the  weathered  material  close  to  the 
bedrock  surface  may  locally  be  permeable.  However,  the  assumption  that 
the  bedrock  is  impermeable  relative  to  the  alluvial  aquifer  is  believed 
to  be  valid.  Evaluation  of  the  possible  minor  flows  within  the  bedrock 
is  not  included  in  the  scope  of  this  study,  nor  does  it  currently  appear 
warranted. 
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A  large  data  base  exists  to  evaluate  the  material  properties  (both  of 
the  bedrock  and  unconsolidated  alluvial  deposits)  within  and  to  some 
degree  in  the  near  vicinity  of  the  Arsenal  grounds.  The  location  of  the 
available  boreholes  and  wells  is  provided  in  Figure  3.  The  same  quality 
and  type  of  information  is  not  available  for  some  of  the  boreholes  and 
wells  shown  in  Figure  3.  Consequently  any  one  of  the  representations 
developed  from  this  data  base  does  not  typically  include  information 
from  all  the  boreholes  or  wells  present  in  Figure  3.  Relative  to 
defining  the  top  of  bedrock  the  majority  of  the  locations  identified  in 
Figure  3  provided  usable  information  that  D'Appolonia  was  able  to 
develop  into  an  accurate  top  of  bedrock  contour  map  (Figure  4) . 

2.2  AQUIFER  PROPERTIES 

For  the  majority  of  the  existing  boring  logs,  the  materials  character¬ 
istics  of  the  unconsolidateddeposits  were  indicated  only  by  a  Unified 
Soil  Classification  System  (USCS)  group  symbol.  The  following  USCS 
groups  were  typically  considered  to  be  aquifers: 

GW  -  well  graded  gravels,  gravel-sand  mixtures, 
little  or  no  fines. 

GP  -  poorly  graded  gravels,  gravel-sand  mixtures, 
little  or  no  fines. 

GM  -  Silty  gravels,  poorly  graded  gravel-sand-silt 
mixtures. 

GC  -  cleyey  gravels,  poorly  graded  gravel-sand-clay 
mixtures. 

SW  -  well  graded  sands,  gravelly  sands,  little  or  no 
fines. 

SP  -  poorly  graded  sands,  gravelly  sands,  little  or 
no  fines 

,SM  —  silty  sands,  poorly  graded  sand— silt  mixtures 

SC  —  clayey  sands,  poorly  graded  sand— clay  mixtures 
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Various  combinations  of  these  groups  such  as 


SPGP  and  SPSM  were  also 


Permeability  values  for  the  sand,  and  sand  and  gravel  units  were  obtained 
from  aquifer  tests  in  1976  (Mitchell.  1976)  and  1978  (Vispi,  1978) . 
D'Appolonia  (1979)  evaluated  these  test  results  and  estimated  representa¬ 
tive  permeabilities  for  the  sand,  and  sand  and  gravel  units  of  about 
3,000  gpd/ft2  (400  ft/day)  and  5,000  gpd/ft2  (668  ft/day),  respect¬ 
ively.  An  order  of  magnitude  reduction  in  permeability  exists  for  units 
that  contain  appreciable  silt.  Any  unit  containing  clay  has  a  negligible 
permeability  compared  with  the  clean  sands,  and  sand  and  gravels. 


The  potential  variance  or  error  present  in  the  USCS  designations  on  the 
Arsenal's  data  is  important  to  note  when  evaluating  boring  logs  that 
have  been  accumulated  over  a  period  as  long  as  20  years.  The  exact 
distinction  of  clayey  or  silty  soils  in  the  USCS  can  frequently  only  be 
made  via  laboratory  tests  (liquid  limit  and  plastic  limit).  In  the 

«:d,  these  parameters  can  be  difficult  to  accurately  quantify  unless 
observer(s)  has  substantial  experience.  Therefore,  the  possibility 
exists  that  a  soil  classified  as  silty  in  one  boring,  for  example,  may 
have  been  described  as  clayey  in  another  boring  by  a  different  inspector. 
In  addition,  relatively  small  amounts  of  clay  or  silt  may  significantly 
affect  permeability  values.  For  these  reasons,  permeability  distribu¬ 
tions  from  the  boring  logs  alone  are  difficult  to  interpret  in  exact 
terms .  , 


2.3  ALLUVIAL  AQUIFER  THICKNESS  AND  DISTRIBUTION 

Figure  5  is  a  plot  of  alluvial  aquifer  thickness  as  indicated  by  the 
existing  boring  logs.  The  aquifer  is  defined  as  all  sand  and  gravel, 
gravel,  and  sand  units  that  are  either  unconfined  or  confined  below  an 
impermeable  layer.  In  both  cases,  the  aquifer  thickness  is  the  thick¬ 
ness  of  permeable  materials  above  and  below  the  water  table.  Permeable 
deposits  lying  above  a  confining  layer  are  included  only  if  partly  or 
wholly  saturated,  the  thickness  counted  being  the  thickness  of  the 
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entire  permeable  zone  above  and  below  the  saturated  zone.  To  evaluate 
the  saturated  thickness,  the  difference  between  the  potent iometric  and 
bedrock  surface  was  used. 

The  thickness  of  the  aquifer  varies  from  zero  on  bedrock  highs  to  a 
recorded  maximum  of  41  feet  about  1  mile  south  of  the  containment 
system.  As  shown  in  Figure  5,  thicknesses  are  generally  greater  in  the 
southern  and  central  portions  of  the  area  of  interest.  The  contoured 
data  illustrates  the  presence  of  two  sediment  troughs  that  correspond 
approximately  with  the  bedrock  valleys  described  above.  In  detail,  many 
locations  show  a  significant  variation  in  thickness  over  a  short  horizon¬ 
tal  distance.  This  variation  is  probably  real  resulting  from  the  lenti¬ 
cular  nature  of  coarser  grained  channel  fill  deposits  within  the  aquifer. 
In  some  cases,  however,  the  variation  may  be  due  in  part  to  differing 
interpretations  of  the  material  properties  of  the  same  deposits  during 
logging.  The  use  of  smoothed  contours  highlights  the  major  trends  in  the 
aquifer  rather  than  minor  channeling  effects. 

Working  cross-sections  were  prepared  at  various  locations  across  the 
aquifer  normal  to  the  direction  of  groundwater  flow  to  assess  aquifer 
continuity.  These  sections  suggest  that  the  aquifer  is  relatively 
continuous.  While  the  thickness  does  vary  locally  there  are  not  signif¬ 
icant-continuous  impermeable  barriers  between  the  deeper  permeable 
channel  fill  deposits  that  could  have  a  significant  effect  on  flow 
direction  and  distribution. 

ft 

2.4  PERMEABILITY  DISTRIBUTION 

As  noted  in  Section  2.2,  silty  units  have  a  significantly  lower  perme¬ 
ability  than  the  clean  sands  or  sands  and  gravels.  Lithologic  maps 
were  prepared  for  working  purposes  to  show  the  predominant  character  of 
the  alluvial  aquifer.  These  maps  show  a  tendency  for  a  higher  propor¬ 
tion  of  silty  units  towards  the  east  and  within  the  deeper  aquifer 
channel  fill  deposits.  In  all  cases,  when  there  is  a  high  proportion 
of  silty  units,  however,  adjacent  borings  show  a  high  proportion  of 
sand  or  sand  and  gravel.  This  observation  suggests  that  there  is  a 
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sufficient  proportion  of  sand  or  sand  and  gravel  throughout  the  entire  aquifer. 
This  result  indicates  that  the  aquifer  can  reasonably  be  modeled  on  a  relative¬ 
ly  macro  scale  to  successfully  predict  the  total  flux  across  the  Arsenal  s 
northern  boundary.  Locally,  a  concentration  of  flow  within  restricted  channels 
of  higher  permeability  would  be  expected  to  be  present.  It  is  not,  however, 
practical  to  model  these  channels  given  the  existing  data  (despite  the  large 
number  of  borings)  and  the  complex  micro  flow  systems. 

2.5  POTENTIOMETRIC  SURFACE 

The  average  spring  prepilot  plant  potentiometric  surface  for  the  spring  1978 
is  provided  in  Figure  6a.  The  potentiometric  gradient  is  generally  toward 
the  north  with  an  average  gradient  of  0.006  ft. /ft.  (Kolmer  and  Anderson, 

1977) .  The  observed  gradient  in  the  bedrock  valley  between  the  north  boundary 
and  Reservoir  "F"  is  much  lower.  This  observation  suggests  that  the  permeability 
^^^hin  this  zone  is  much  higher  than  the  material  along  the  north  boundary 
^Plre  the  gradient  is  steeper. 

Alluvium  water  levels  from  May  to  June  1979  (Figure  6b)  gathered  by  Zabell 
(1979)  demonstrate  similar  flow  characteristics.  Gradients  in  the  range  of 
0.006  to  0.008  ft/ft.  along  the  north  boundary  are  evident.  Comparison  of  the 
data  reveal  significant  similarities  suggesting  the  alluvial  groundwater  system 
is  presently  at  or  near  steady  state  conditions. 


13 


lO^^IPCDILCDFim 


Alluvium  Water  Levels  (May-June  1979) 
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■ASSESSMENT  A: 


D'Appolpnia  Consulting  Engineers 


To  develop  an  optimum  design  for  the  dewatering  and  recharge  subsystems 
of  the  north  boundary  containment  system,  a  conceptual  hydrogeologic 
model  was  formulated  into  a  finite  element  mathematical  model.  The 
objective  of  this  modeling  effort  was  to  evaluate  natural  subsurface 
outflow  across  the  RMA's  north  boundary  and  to  develop  performance 
predictions  for  the  designed  system.  The  combination  of  a  vari¬ 
able  potent iometric  gradient,  contrasts  in  aquifer  permeability 
and  the  irregular  boundary  conditions  suggest  that  the  use  of  a  mathe¬ 
matical  flow  model  is  desirable  to  reflect  the  detailed  nature  of  the 
available  data. 


3.1  MODEL  DESCRIPTION 

The  computer  program  used  in  this  study  is  a  two-dimensional  finite 
element  code  which  simulates  the  performance  of  an  aquifer  on  a  regional 
basis.  The  basic  governing  equations  of  the  flow  regime  are  well  estab- 
^jj^ed  (Bear,  1972;  Finder  and  Frind,  1972),  and  described  later  in  this 
section.  To  model  the  performance  of  the  north  boundary  alluvial  aquifer, 
the  region  was  divided  into  a  system  of  grids  which  are  called  elements. 
When  hydrodynamic  parameters  and  stresses  are  provided  (such  as  trans¬ 
missivity,  storage  coefficients,  pumping  and  recharge  rates),  the  program 
calculates  time-variable  or  steady  state  potentiometric  surfaces  and 
consequently  the  flow  vectors  of  the  region.  A  full  description  of  the 
methodology  used*,  assumptions  made,  and  accuracy  in  the  program  is 
contained  in  Appendix  A.  A  summary  of  this  information  follows. 

Governing  Equations 

The  combined  equation  of  motion  and  continuity  for  flow  in  a  two-dimen¬ 
sional  horizontal  plane  can  be  written; 


>r  3<l>  j.  J 7  9<l>  _  p  =  c  .li 

ax  XX  ax  ^  ay  ^yy  ay  ^  at 
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where : 

''^xx  “  principal  hydraulic  transmissivity  along  x-direction 

[L2/t]. 

Tyy  =  principal  hydraulic  transmissivity ’ along  y-direction 
,[L2/t]. 

♦  =  potentiometric  surface  [L]. 

P  =  the  strength  of  a  sink  or  source  [L3/L2t]. 

S  =  storage  coefficient  [dimensionless]. 
x,y  =  Cartesian  coordinates  [L] . 
t  =  time  [t] . 


8  a  a 

lx  Ty  "ait  ~  partial  derivatives  with  respect  to  x,  y,  and 
’  t,  respectively. 


The  transmissivity  of  an  aquifer  under  artesian  conditions  is  the  product 
of  the  hydraulic  conductivity  (permeability)  and  the  thickness  of  that 
aquifer.  Thus, 


"^xx  “  ^  ^xx>  "^yy  “  b  K 


where : 


b  =  thickness  of  aquifer  [L] . 


^yy 


(2) 


Kjjx  =  principal  hydraulic  conductivity  along  x-direction 
[L/t]. 


Kyy  =  principal  hydraulic  conductivity  along  y-direction 
.  [L/t]. 


Directional  anisotropy  in  permeability  was  not  implemented  in  the 
modeling  of  the  north  boundary  alluvial  aquifer.  In  a  water  table 
aquifer,  the  saturated  thickness  is  used  to  calculate  the  transmissi¬ 
vity  in  an  iterative  fashion.  This  process  has  the  effect  of  decreas¬ 
ing  transmissivity  in  proportion  to  drawdown.  A  more  complete  descrip¬ 
tion  of  the  governing  equations  is  provided  in  Appendix  A,  Section  1.4. 
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Basic  Asstimptions 

^he  following  assumptions  within  the  model  are  valid  for  regional  ground- 
water  flow: 


a.  The  flow  is  essentially  horizontal  in  a  twor 
dimensional  plane.  This  assumption  is  valid  when 
the  variation  of  thickness  of  the  aquifer  is  much 
smaller  than  the  thickness  itself.  This  approxi¬ 
mation  fails  in  regions  where  the  flow  has  a 
vertical  component. 

b.  The  fluid  is  homogeneous  and  slightly  compressible. 

c.  The  aquifer  is  elastic  and  generally  nonhomogen- 
eous  and  anisotropic.  The  consolidating  plastic 
medium  deforms  during  flow  due  to  changes  in 
effective  stress  with  only  vertical  compress¬ 
ibility  being  considered. 

d.  For  the  two-dimensional  horizontal  flow  assumption, 
an  integrated  potentiometric  level  is  used  where 
the  value  is  determined  along  vertical  lines 
extending  from  the  bottom  to  the  top  of  the 
aquifer. 

e.  Within  each  element,  parameters  such  as  trans¬ 
missivity,  and  storage  coeficient  are  assumed  to 
remain  constant  within  each  time  step. 


Method  of  Solution 

The  Galerkin  finite  element  technique  is  used  to  simultaneously  solve 
for  flow  and  potentiometric  head  within  the  model.  The  direct  method  of 
solution  of  the  simultaneous  equations  of  flow  are  used  to  avoid  errors 
associated  with  iterative  methods.  A  detailed  description  of  the 
method  of  solution  is  provided  in  Appendix  A,  Section  1.5. 

Initial  Boundary  Conditions 

To  solve  the  equation,  certain  additional  conditions  imposed  by  the 
physical  situation  at  the  boundaries  are  selected.  In  general,  for  flow 
through  an  aquifer,  three  different  boundary  conditions  are  applicable: 
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1.  Known  potent iometric  level  on  the  boundary; 

2.  Known  flux  on  the  boundary;  and 

3.  Known  potent iometric  level  and  its  normal  deriva¬ 
tive  on  the  boundary. 

Any  of  the  above  boundary  conditions  or  their  combinations  may  appear  in 
the  modeling. 

At  the  initial  time,  either  the  potent iometric  surfaces  are  known  in 
the  entire  domain  or  the  hydrologic  stresses  (such  as  pumping  and 
recharge)  are  specified  and  boundary  conditions  are  known.  For  the 
second  case,  the  system  has  reached  the  steady-state  condition;  there¬ 
fore,  the  initial  potent iometric  surfaces  can  be  found  by  solving 
Equation  (1),  while  setting  Ml^t  equal  to  zero.  A  more  complete 
description  of  the  initial  and  boundary  conditions  may  be  found  in 
Appendix  A,  Section  1.4.3.  The  accuracy,  limitations  and  restrictions 
of  the  model  are  detailed  in  Appendix  A,  Section  1.6. 

3,2  MODEL  INPUT  DATA  REQUIEEMENTS 

To  solve  Equation  (1),  the  hydrodynamic  parameters  and  stresses  must  be 
specified,  and  initial  and  boundary  conditions  assigned.  The  required 
input  data  to  solve  the  regional  groundwater  flow  in  an  aquifer 
follow: 

o  Type  of  Aquifer 
o  Static  Potentiometric  Surface 
o  Pumping  Rate 
o  Hydraulic  Conductivity 
o  Thickness  of  Aquifer 
o  Storage  Coefficient 
o  Boundary  Conditions 

Along  with  the  above  information,  any  expected  alterations  in  the 
groundwater  flow  regime  (such  as  lowering  the  potentiometric  surface) 
should  also  be  incorporated  into  the  model. 
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3.3  ■  MODELING  PROCEDURE 
Grid  System  Development 

9  egion  around  the  RMA's  north  boundary  was  divided  into  grids  of 
quadrilateral  and  triangular  elements  as  shown  in  Figure  7.  This  grid 
system  consists  of  392  elements  and  437  nodes.  To  increase  the  utility 
of  the  model  as  a  design  tool,  the  grid  mesh  is  finer  in  the'  area 
of  the  north  boundary  containment  system".  The  grid  size  was  expanded 
in  areas  of  sparse  data  and  near  the  boundaries.  Input  data  were 
prepared  predominantly  in  the  form  of  overlay  maps.  Discrete  values 
for  the  various  parameters  were  coded  into  the  model  at  each  element. 


Type  of  Aquifer 

As  discussed  in  Chapter  2.0,  the  aquifer  of  concern  in  the  north 
boundary  area  consists  of  a  shallow,  unconsolidated  alluvial  aquifer. 
Depending  on  the  location,  this  aquifer  acts  as  a  phreatic  or  water 
table  aquifer  with  locally  confined  areas.  The  areas  where  the  poten- 
tiometric  surface  rises  above  the  top  of  the  main  productive  unit 
typically  do  not  react  strictly  as  a  confined  aquifer  due  to  the  silty 
na^re  of  the  overlying  units.  The  typically  lenticular  nature  of  the 
a^jj^ter  causes  this  aquifer  to  respond  as  a  semi-confined  aquifer 
during  short-term  stress  and  as  a  water  table  aquifer  after  a  longer 
period  of  time.  Recharge  to  the  aquifer  is  primarily  due  to  leakage 
from  surface  water  bodies.  Areal  recharge  due  to  infiltration  of 
precipitation  is  assumed  to  be  negligible  due  to  the  semi-arid  nature 
of  the  climate.  The  lower  boundary  of  the  alluvial  aquifer  is  assumed 
to  be  impervious  relative  to  the  alluvium. 


Boundary  Conditions 

Two  types  of  boundary  conditions  were  assumed  around  the  periphery  of  the 
finite  element  grid  system,  namely  either  a  constant  head  boundary  or  a 
no  flow  boundary.  Figure  7  illustrates  the  location  where  both  of  these 
boundary  types  were  used.  Constant  head  boundaries  were  imposed  at  the 
upgradient  ends  of  the  alluvium  filled  paleovalleys  that  feed  the  north 
boundary  alluvial  aquifer.  Additional  constant  head  areas  were  set  at 
the  downgradient  end  of  the  grid  system  to  simulate  regional  outflow. 
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Physical  boundaries  of  the  alluvial  aquifer,  generally  associated  with 
pinchouts  over  bedrock  highs  reported  by  Robson  (1976),  were  treated  as 
no  flow  boundaries.  The  net  effect  of  this  set  of  boundary  assumptions 
is  to  simulate  regional  inflow  from  the  south  and  outflow  at  the  north 
end  of  the  grid  system. 

Hydrodynamic  Parameters 

The  hydrodynamic  parameters  of  the  aquifer  system  utilized  in  the  final 
computer  runs  presented  in  this  section.' |  are  as  follows: 

Unit  Permeability  Storage  Coefficient 

_ (ft. /day)  _ 

Impermeable  areas  0.1  .10 

Silty  sands  (SM)  53.47  .10 

Sands  and  gravel  (SPGP)  401.05  .10 

Slightly  silty  sands  (SMSP)  200.52  .10 

These  permeabilities  were  based  on  evaluations  of  the  results  of  the 
two  field  aquifer  testing  programs  cited  in  Section '2.0,  interpretation 
from  boring  logs,  analysis  of  grain  size  curves,  and  the  result  of 
interactive  model  calibration. 

The  effective  thickness  of  the  aquifer  is  based  on  the  calculated 
difference  between  the  bedrock  elevation  and  the  potentiometric  surface 
within  an  element  and  is  updated  during  each  time  step.  The  values  for 
saturated  thickness  were  checked  by  the  program  so  that  saturated 
thickness  does  not  exceed  the  thickness  of  the  potentially  productive 
aquifer.  The  values  used  for  bedrock  elevation  and  aquifer  thickness 
were  presented  in  Section  2.0. 


3.4  COMPUTATION  METHODOLOGY 

The  initial  model  calibration  and  computation  of  steady  state  potentio¬ 
metric  levels  were  conducted  to  assess  the  accuracy  of  the  model. 
Initial  estimates  of  permeability  within  each  of  the  zones  identified 
in  Figure  5  were  based  on  the  field  testing  data  available.  The 
constant  head  boundaries  were  set  at  the  observed  values  from  the 
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pote'ntiometric  surface  maps  shown  in  Figure  6.  The  model  was  then 
exerted  to  obtain  steady  state  conditions  with  the  resulting  calcu- 
laWpotentiometric  surface  being  compared  to  the  observed  values. 
Permeabilities  were  adjusted  and  the  steady  state  potent iometric  levels 
were  recalculated.  This  calibration  proceeded  in  an  iterative  fashion 
until  the  agreement  between  calculated  and  observed  levels  was  judged 
to  be  adequate.  The  final  steady  state  calculated  potentometric  level 
map  is  shown  in  Figure  8.  The  general  agreement  between  calculated  and 
observed  potentiometric  levels  is  adequate  for  purposes  of  evaluation 
of  well  field  design  and  total  flux  in  the  aquifer  across  the  RMA's 
north  boundary.  Some  of  the  details  of  the  potentiometric  surfa«  are 
not  duplicated,  particularly  in  the  southwest  quadrant  of  the  grid 
network  where  gradients  are  low.  The  steady  state  levels  shown  in 
Figure  8  were  used  as  the  initial  conditions  in  subsequent  model  runs; 
therefore,  the  groundwater  system  began  in  a  state  of  equilibrium  for 
design  evaluation. 


Steady  state  flows  across  the  containment  system  area  were  evaluated 

u^^  the  computer  model.  The  y  component  of  flow  was  calculated  for 

eW  element  across  the  grid  system  (element  numbers  178  through  198, 

shown  in  Figure  7  using  the  relationship  q  =  31.  This  equation  was 

ay 

numerically  evaluated  at  the "center  of  each  element  and  when  multiplied 
by  the  element  width,  the  total  flow  was  estimated.  The  flow  rates  by 
element  are  summarized  in  Table  1. 

Simulation  of  the  containment  system  vas  based  on  the  proven  design 
concept  demonstrated  during  operation  of  the  pilot  containment  system. 
This  design  consists  of  the  following. 

o  A  slurry  wall  to  block  the  flow  of  groundwater; 

o  Dewatering  wells  located  in  the  upgradient 
direction,  and 

o  Injection  wells  located  in  the  downgradient 
direction. 
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The  objective  of  this  system  is  to  operate  with  as  little  overall  impact 
on  water  levels  as  possible.  This  condition  is  accomplished  by  operating 
the  well  field  so  that  the  quantity  of  water  pumped  is  equal  to  that  which 
would  normally  flow  across  the  boundary  if  no  barrier  were  present. 

The  presence  and  operation  of  this  system  was  simulated  by  changing  the 
permeability  of  one  row  of  elements  to  a  very  low  value  along  a  line  cor¬ 
responding  to  the  current  alignment  of  the  pilot  system  wall.  Wells  were 
simulated  300  feet  upgradient  from  the  wall  with  pumping  rates  equal  to  the 
computed  steady  state  flux  within  each  element.  A  similar  setup  was  established 
downgradient  from  the  wall  for  the  recharge  wells.  Operation  of  this  system 
was  simulated  for  a  period  of  300  days  using  a  time  step  of  10  days  for 
accuracy.  The  resulting  average  potentiometric  surface  configuration  is 
shown  in  Figure  9.  The  levels  represented  in  this  figure  are  averages 
within  each  element.  The  levels  within  the  dewatering  wells  themselves  will 
be  somewhat  lower.  In  the  recharge  wells,  the  levels  will  be  higher  than 
shown.  At  distances  of  about  150  feet  from  the  wells,  however.  Figure  9 
is  representative  of  anticipated  conditions. 

3.5  TOTAL  QUANTITIES  OF  CONTAMINANTS  TO  BE  TREATED 

A  contoured  concentration  map  of  fluoride  has  been  generated.  The  data  for 
this  map  was  provided  by  RMA  personnel.  This  F  map  was  prepared  using  data 
collected  during  May  to  August  1979  and  contained  in  the  Material  Analysis 
Laboratory  Report  submitted  to  the  Geohydrology  Division,  September  21, 

1979.  Concentrations  on  the  map  represent  the  results  for  samples  collected 
after  the  pilot  plant  began  operation  on  July  28,  1978.  The  concentration 
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data  for  fluoride  was  plotted  at  the  appropriate  borehole  locations  (Figure  3) . 
This  information  was  then  contoured  using  linear  interpolation  between  known 
data  points.  Variance  to  the  linear  iiiterpolation  procedure  occurred  only 
in  a  few  localized  instances  where  a  temporal  comparison  with  other  maps 
of  the  same  contaminant  suggested  a  more  realistic  contour  location.  Ground- 
water  dispersion  of  chemical  species  is  not  expected  to  be  a  linear  function. 
Therefore,  contours  located  in  areas  of  sparce  data  may  not  represent  actual 
concentrations  in  that  area.  However,  in  most  cases,  the  amount  of  data  available 
is  enough  to  insure  reliable  conclusions. 

The  resulting  contoured  map  (Figure  lo)  reveals  that  there  are  three  high 
concentration  areas.  One  high  is  centered  at  the  base  of  the  western  bedrock 
^^h.  The  other  two  highs,  are  located  slightly  to  the  east.  The  concentration 
^^tours  have  a  west  to  east  trend  Chigh  to  low) .  Over  60  percent  of  Section 
24  has  a  concentration  of  <2  mg/il.  Relatively  low  concentration  areas  exist 
in  areas  of  high  DCPD  and  DBCP  concentrations.  This  result  probably  indicates 
that  much  of  the  fluoride  has  been  removed  by  the  large  flow  of  groundwater 
through  these  areas.  The  operation  of  the  pilot  containment  system  has 
had  little  effect  upon  the  location  of  fluoride  concentration  bands. 


Figure  10  was  generated  using  data  collected  during  May  to  August  1979. 

Maps  were  also  generated  using  data  collected  during  the  winter  months 
(October  and  December  1978).  These  maps  were  evaluated  for  comparative 
purposes  and  are  not  included  in  this  report.  The  general  trends  and  locations 
of  the  contours  were  again  similar.  In  most  cases,  the  seasonal  variations 


fluoride  concentrations  were  minimal.  However,  certain  locations  reveal 
large  differences:  Borehole  168,  winter  =  .68  mg/ 1  fluoride,  summer  =  4.6 
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mg/il  fluoride;  and  Borehole  539,  winter  =  0.36  mg/i  fluoride,  summer  = 

1.5  mg/Jl  fluoride. 

Calculation  of  Quantities  of  F~  to  be  Processed 

The  concentration  in  the  center  of  Elements  178-198  was  estimated  from 
the  contoured  concentration  maps  (Figure  10)  for  F  .  Flow  rates  | 

used  in  calculational  loadings  were  derived  from  the  steady  state  model 
results  described  in  Section/ 3.0.  The  center  of  these  elements  are  located 
approximately  200  feet  upgradient  from  the  location  of  the  dewatering  wells 
in  the  pilot  system.  Therefore,  the  recovery  values  calculated  will  pre¬ 
dict  the  recovery  after  operation  has  started  and  equilibrium  is  reached. 

The  pumping  rate  for  each  element  is  also  known;  therefore,  the  total 
weight  of  F“  contained  in  the  water  can  be  calculated.  These  values  are 
provided  in  Table  2:/  Approximately  6.6  kg  of  fluoride  ion  (F")  will  be 
handled  each  day.  This  corresponds  to  an  expected  average  concentration 
in  all  elements  of  2.72  mg/i.  Because  of  the  large  amount  of  data  available, 
the  weights  of  F“  processed  should  be  fairly  accurate. 

Changes  in  the  Quantities 

Contour  maps  of  the  total  area  of  Sections  23  and  24  were  generated  so 
that  any  future  alternation  in  the  concentrations  of  the  contaminants 
could  be  predicted.  Basically,  the  contour  lines  run  parallel  to  the 
water  flow.  Therefore,  concentrations  will  vary  little  with  time. 

However,  there  will  be  an  overall  decrease  in  the  concentrations  as  the 
aquifer  is  flushed.  Where  contour  intervals  run  east  and  west,  the 
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upgradient  value  is  generally  lower,  therefore  the  overall  concentrations 
will  probably  decrease. 
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TABLE  1 

STEADY  STATE  FLOWS  ACROSS  THE  NORTH  BOUNDARY 


Element  Number^ 

Flow  Rate(2) 

178 

3 

179 

14740 

180 

32105 

181 

36287 

182 

27400 

183 

23043 

184 

42347 

185 

33542 

186 

42354 

187 

39676 

188 

41045 

189 

45780 

190 

39856 

191 

34279 

192 

30774 

193 

66246 

194 

63284 

195 

25632 

196 

0 

*  197 

0 

198 

0 

TOTAL  FLOW  RATE  638393 

NOTES: 

(l^See  Figure  7  for  locations  of  elements. 
(2)piow  rates  reported  in  gallons  per  day  (gpd) 
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TABLE  2 


TOTAL  QUANTITIES  OF  FLUORIDE  TO  BE  PROCESSED 


CONCENTRATION 

GRAMS /DAY 

FLOW  .RATE 

- fT- - 

(GPD) 

(mg/A) 

F- 

3 

4.4 

0 

14740 

5.4 

301 

32105 

5.6 

681 

36287 

3.8 

522 

27400 

4.0 

415 

23043 

4.5 

393 

42347 

4.7 

753 

33542  ■ 

4.1 

521 

42354 

3.4 

545 

39676 

2.9 

436 

41045 

2.6 

404 

45780 

2.4 

416 

39856 

2.1 

317 

34279 

1.8 

234 

30774 

1.4 

163 

66246 

1.0 

251 

63284 

0.7 

168 

25632 

0.5 

49 

0 

0.5 

0 

0 

0.7 

0 

0 

1.0 

0 

ELEMENT 
NO. 


178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 
1 

92 

193 

194 

195 

196 

197 

198 


m 


TOTAL 


638393 


6566 
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4.0  ASSESSMENT  B:  US  Army  Toxic  and  Hazardous 

Materials  Agency 

To  determine  the  expected  contaminant  loading  of  a  compound  within  a 
controlled  section  of  aquifer,  a  mass  flux  computation  was  formulated. 

By  knowing  total  groundwater  flow  and  total  contaminant  mass  migrating 
per  unit  time  in  any  section  of  aquifer,  a  concentration  prediction 

can  be  made.  This  technique  addresses  variable  aquifer  characteristics  / 
such  as  cross  sectional  area  (saturated  thickness  and  length) , 
potentiometric  gradient  and  aquifer  permeability. 


4.1  MASS  FLUX  DESCRIPTION 

The  computer  program  used  in  this  study  utilizes  a  two-dimensional 
jjj^resentation  of  an  aquifer  to  model  groundwater  flow  characteristics. 
The  basic  equation  uses  a  form  of  Darcy’s  law  for  the  flow  of  water 
through  a  porous  media.  To  compute  the  contaminant  mass  flux  for  the 
north  boundary  alluvial  aquifer,  a  transect  perpendicular  to  the 
northward  flowing  groundwater  was  established  and  divided  into  a 
series  of  rectangular  elements.  Each  element  was  then  described 
with  appropriate  geohydrologic  terms  to  thoroughly  define  its  flow 
characteristics.  Past  contaminant  concentration  data  was  assessed 


to  determine  a  representative  concentration  value  for  that  section 
of  aquifer.  Input  of  this  information  for  every  element  yielded  a 
composite  mass  flux  estimate.  A  full  description  of  the  software 
input /output  requirements  is  contained  in  Appendix  B.  A  summary  of 


;he  governing  equation,  assumptions  made  and  data  used  in  the 
ogram  follows . 
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Governing  Equation 


The  mass  flux  computation  uses  a  combination  of  groundwater  flow  quantity 
(volume  per  time)  and  contaminant  concentration  profile  (weight  per  volume) 
to  calculate  contaminant  mass  movement  (weight  per  time).  A  form  of  Darcy's 
law  for  the  flow  of  groundwater  forms  the  basis  of  the  mass  flux  flow  estimate. 
The  governing  equation  thus  can  be  represented  by: 

0  -  KLtI  X  C  X  U  (3) 

where : 

0  =  mass  flux  across  a  transect  in  gms/day 
K  =  permeability  of  saturated  zone  in  feet/day 
L  =  saturated  length  in  feet 
t  =  saturated  flow  thickness  in  feet 
I  =  gradient  of  the  water  table 

C  =  contaminant  concentration  in  mg/liter  or  ug/ liter 
TJ  =  units'  adjustment 


Basic  Asstimptions 

The  following  assumptions  are  necessary: 

a.  The  flow  is  essentially  horizontal  in  a  two-dimensional  plane. 
This  assumption  is  valid  when  the  variation  of  thickness  of 

of  the  aquifer  is  much  smaller  than  the  thickness  itself. 

This  approximation  fails  in  regions  where  the  flow  has  a 
vertical  component. 

b.  The  fluid  is  homogeneous  and  slightly  compressible. 

c.  For  the  two-dimensional  horizontal  flow  assumption, 
an  integrated  potentiometric  level  is  used  where 

the  value  is  determined  along  vertical  lines  extending 
from  the  bottom  to  the  top  of  the  aquifer. 
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d.  Within  each  element,  parameters  such  as  cross-sectional 
areas,  permeability  and  potentiometric  gradients  are 
assumed  to  remain  constant.  This  assumption  dictates 
that  the  aquifer  under  consideration  is  at  steady  state. 


4.2  MASS  FLUX  INPUT  DATA  REQUIREMENTS 

To  solve  Equation  (3)  ,  the  groundwater  flow  characteristics  must  be  specified 

for  a  section  of  aquifer  and  initial  boundaries  assigned  for  the  contaminant 

concentration  assessment.  The  required  input  data  to  solve  the  mass  flux 

computation  follow: 

o  Transect  Specifications 
o  Potentiometric  Gradients 
o  Aquifer  Cross-sectional  Area  (Saturated) 
o  Aquifer  Permeability 
o  Contaminant  Concentration  Assessment 

l^^ual  input  values  with  rationale  for  each  of  these  areas  is  described  in 

the  following  subsections. 


Transect  Specifications 

Because  of  the  assumption  of  two-dimensional  flow,  the  transect  through  which 
the  mass  flux  is  to  be  computed  must  be  specified  perpendicular  to  the  pre¬ 
dominant  groundwater  flow  direction.  Intervals  along  the  transect  are  chosen 
to  define  sides  of  the  various  rectangular  elements  in  which  groundwater  flow 
will  be  calculated.  The  larger  the  number  of  intervals  along  a  transect,  the 
greater  the  mass  flux  refinement  can  be.  However,  directly  associated  with 
this  increase  in  intervals  is  the  enlargement  of  geohydrologic  data  needed  to 
define  each  aquifer  element. 


For  the  situation  of  the  north  boundary  of  BMA.,  a  two  mile  transect  parallel  ] 
to  the  north  boundary  was  chosen  for  analysis  (Figure  11) .  This  transect  lies 
approximately  500'  back  from  the  arsenal  boundary  (in  line  with  the  pilot  system) 
in  the  north  part  of  sections  23  and  24.  The  transect  was  subsequently 
divided  into  250'  intervals  which  resulted  in  approximately  40  equally  spaced  \ 
aquifer  units.  Groundwater  flow  in  this  area  is  at  perpendicular  direction 
to  the  boundary. 

Potentiometric  Gradients 

Aquifer  gradients  must  be  established  for  each  transect  interval.  Potentio¬ 
metric  gradients  represent  the  driving  force  for  the  groundwater  system. 

Historical  water  level  data  was  reviewed  to  determine  possible  long-term 
trends  in  a  alluvial  aquifer.  New  borings  and  piezometers  were  located  to 
fill  gaps  in  the  existing  boring  and  piezometer  array.  Data  from  the  logs 
and  piezometers  were  used  to  determine  water  table  elevations,  base  of  alluvial 
aquifer  and  groundwater  flow  patterns.  A  general  overview  of  this  data 
has  been  previously  highlighted  in  Section  2;  Site  Hydrogeology.  Specific 
hydrodynamic  parameters  for  input  in  the  mass  flux  computation  for  the  north 

boundary  area  is  presented  below: 

a.  Water  table  measurements  were  taken  in  the  study  area 
in  October— November  1978  and,  upon  completion  of  the 
installation  of  new  piezometers,  arsenal-wide  measurements 
were  taken  in  March-April  1979  and  in  May- June  1979.  The 
three  sets  of  water  table  measurements  were  reviewed. 

Water  table  fluctuations  during  the  relatively  short 
period  of  this  May- June  1979,  was  constructed  (Figure  6b). 
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The  water  table  contours  show  three  spacing  and  direction 
patterns.  The  first  pattern  occurs  in  the  western  part  of  the 
study  area  (section  23)  and  is  characterized  by  widely  spaced 
contour  lines  (low  gradient)  that  trend  to  the  northeast. 

e  second  pattern  occurs  in  the  eastern  part  of  the  studv  J 
area  (sections  24  and  19).  This  pattern  is  characterized^by 
closely  spaced  contours  trending  in  a  northwest  direction. 
About  1200  feet  south  of  the  north  boundary,  the  first  and 
second  patterns  converge  to  form  closely  spaced  contour 
patterns  that  trend  in  a  northemly  direction. 


c.  Potentiometric  gradients  in  this  area  tend  to  the  north 
with  an  average  gradient  of  0.007  ft/ft.  Small  variations 
in  the  gradient  exist  along  the  boundary  as  one  proceeds  from 
west  to  east.  Estimates  of  this  variability  about  0.007 
have  been  considered  in  the  input  parameters  to  the  mass 
flux  computation  (Table  3). 


Aquifer  Cross-Sectional  Area  (Saturated) 

Ration  3  computes  aquifer  cross-sectional  area  by  the  multiplication  of 
saturated  length  (1)  times  saturated  flow  thickness  (t) .  This  information 


must  be  specified  for  each  transect  interval. 


Borings  placed  in  the  north  boundary  vicinity  have  gathered  geotechnical 
data  to  define  aquifer  media  characteristics.  Previous  boring  locations  and 
logs  were  reviewed  and  a  boring  and  sampling  program  was  designed  to  fill 
data  gaps  in  sections  23,24,25  and  26  and  to  evaluate  potential  pollution  migration 
in  sections  19,23,24,25  and  26.  Seventy-five  new  borings  (Nos.  900-974)  and 
borings  378-380,  382,  385,  and  533  were  located  to  supplement  the  existing 
boring  data.  Split-spoon  samples  were  obtained  at  5  feet  intervals  and  at 
stratum  changes,  where  possible,  from  each  boring.  Samples  were  field  classified 
site  by  several  inspectors.  Laboratory  classification 


was  performed  on 


selected  soil  samples.  Monitorine  wells  were  installed  at  all  new  boring  lo¬ 
cations  for  water  level  measurements,  water  sampling,  and  permeability  tests. 

All  available  logs  and  water  depth  readings  from  piezometers  were  used  to 
construct  the  base  of  the  alluvial  aquifer  CFigure  4) ,  saturated  thickness 
(Figure  5)  and  water  table  (Figure  6a  and  b)  maps  presented  in  section  2.0. 

The  cross  section  at  the  north  boundary  is  shown  in  Figure  12 'On  the  plate, 
two  types  of  data  are  presented— one  cross  section  contains  general  soil  types 
and  water  levels  and  below  it  another  cross  section  presents  the  stratigraphy. 
Specific  hydrodynamic  parameters  for  input  in  the  mass  flux  computation  for 

the  north  boundary  area  is  presented  below: 

a.  Base  of  the  alluvial  aquifer  is  presented  in  Figure  4. 

The  map  generally  depicts  the  weathered  shale  surface 

of  the  Denver  formation  which  underlies  the  alluvium  but 
includes  the  Denver  formation  sands  where  they  are  in 
direct  contact  with  the  alluvium.  This  lower  boundary  is 
assumed  to  be  relatively  impervious.^  The  general  slope 
of  the  base  of  the  alluvial  aquifer  is  to  the  north- 
northeast  north  of  Basin  F  and  north-northwest  in 
section  24.  These  slopes  dictate  the  flow  of  alluvial 
groundwater . 

b.  Figure  5  \is  an  isopach  map  showing  the  saturated  thick¬ 
ness  of  the  alluvial  aquifer  which  includes  Denver  sands 
that  are  in  contact  with  the  alluvium.  The  map  reflects 
saturated  thicknesses  based  on  the  differences  between 
the  base  of  the  alluvial  aquifer  surface  and  piezometer 
surface.  Figure  5  presents  all  saturated  sediments 
which  includes  fine-grained  materials  of  relatively 

low  permeability.  As  expected,  areas  of  greatest  saturated 
thicknesses  follow  channels  and  the  areas  of  least 
saturated  thicknesses  generally  coincide  with  the 
Denver  formation  highs. 

c.  The  mass  flux  transect  of  interest  (cross  section  repre¬ 
sented  in  Figure  12)  is  just  south  of,  and  parallel  to, 
the  north  boundary  and  cross  'sections  23  and  24;  the 
section  is  perpendicular  to  the  northernly  flow  of  ground- 
water  in  the  alluvial  aquifer.  This  section  identifies  the 
base  of  the  alluvial  aquifer  as  weathered  Denver  shale 
except  for  two  areas  which  are  identified  later.  Denver 
High  "A"  occurs  near  the  center  of  the  western  half  of  the 
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cross  section  and  is  the  northemly  extension  of  the  Denver 
High  "A"  in  cross  section  B-B*.  West  of  the  Denver  High  "A" 
the  Denver  surface  decreases  in  elevation  and  contains  two 
small  channels.  Between  these  two  channels,  the  Denver 
formation  consists  of  silt  and  clayey  sand.  East  of  the 
High  "A"  is  one  small,  narrow  channel  and  one  shallow, 
wide  channel  followed  by  a  wide,  flat  surface  extending 
across  the  boundary  of  sections  23  and  24.  The  center  of 
section  24  has  one  large  channel  (identified  in  Figure  12 
as  channel  1)  with  two  smaller  channels  to  its  west.  The 
Denver  formation  between  the  channels  consists  of  silty 
and  clayey  sands.  East  of  channel  1  the  Denver  surface 
rises  rapidly  and  forms  the  Denver  High  "D".  The  Verdos 
sands  and  gravels,  west  of  the  Denver  High  "A,"  are  from 
0  to  10  feet  thick  and  thicken  to  the  west.  East  of  the 
Denver  High  "A"  the  Verdos  alluvium  reappears  and  ranges 
in  thickness  from  5  to  20  feet  with  thickest  deposits 
occurring  in  the  channels .  With  the  exception  of  the 
Denver  High  "A,"  coarse  to  fine  gravels  appear  inter¬ 
mittently  throughout  the  alluvium.  Overlying  the  Verdos 
sands  and  gravels  are  clays,  fine  sands,  and  silts  ranging 
in  thickness  from  5  to  20  feet.  These  sediments  are  in 
direct  contact  with  the  Denver  formation  where  the  Verdos 
is  absent  on  the  High  "A."  These  soils  consist  of  eolian 
sands  and  alluvium  transported  from  higher  elevations , 
except  for  the  areas  in  and  adjacent  to  the  First  Creek 
valley  (center  of  section  24  in  Figure  12)  and  the  small  ' 
valley  west  of  the  section  23  and  24  boundary  where  Piney 
Creek  alluvial  clays,  silts,  and  sands  occur. 

From  the  preceeding  information  it  has  been  determined 
that  the  north  boundary  alluvium  is  fully  saturated 
from  approximately  the  midpoint  of  interval  11  to  just 
inside  interval  35.  Saturated  thickness  of  this  aquifer 
ranges  from  1  to  22  feet.  Exact  cross  sectional  data 
used  for  input  are  delineated  in  Table  3. 


Aquifer  Permeability 

Field  pump  tests  and  rising  and  falling  head  (slug)  tests  were  performed 
in  the  study  area  to  determine  the  coefficients  of  permeability  of  the 
alluvial  aquifer.  Specific  hydrodynamic  parameters  for  input  in  the 
mass  flux  computation  for  the  north  boundary  area  is  presented  below. 
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a.  Field  Pxomp  Tests:  Five  pump  tests  were  performed  in  1978 
by  WES.  Three  tests  were  preformed  north  and  northeast 
of  Basin  F  in  section  23  and  two  tests  were  performed 
in  section  24.  One  test  .was  southwest  of  the  sewage 
lagoon  and  the  other  test  was  north  of  the  sewage  lagoon 
near  the  north  boundary.  Wells  345,  368,  529,  548,  and 
549  were  used  for  the  tests.  Observation  wells  were 
installed  on  lines  originating  a±  the  test  well  extending 
to  1000  feet  away  from  the  test  well.  During  the  pump 
tests  water  level  changes  in  the  observation  wells  were 
measured  for  drawdown  and  recovery  and  coefficients 
of  permeability  were  computated  using  the  drawdown  and 
recovery  rates.  Coefficients  of  permeability  computed 
from  test  wells  345,  368,  and  529  ranged  from  2400  to 
12,000  gpd/ft^  and  the  coefficients  of  permeability 
on  the  observation  well  lines  ranged  from  3400  to  8200  gpd/ft'’ 
These  wells  are  in,  or  adjacent  to,  a  subsurface  channel 
which  runs  northeast  from  north  of  Basin  F  towards  the 
north  boundary.  Well  548,  located  astride  a  small  ridge¬ 
like  area,  had  a  coefficient  of  permeability  of  1100 
gpd/ft^;  coefficients  of  permeability  on  the  observation 
well  lines  were  from  1300  to  2000  gpd/ft^.  Well  549, 
located  adjacent  to  Denver  High  "B"  and  in  an  area  of 
rapid  groundwater  gradient  changes,  has  a  coefficient  of 
permeability  of  430  gpd/ft^.  The  observation  well 
lines  reflected  rapid  changes  in  the  coefficient  of 
permeability  when  measured  over  a  short  distance 
(0-50  ft;  250  gpd/ft^)  as  compared  with  a  long  distance 
(50-1000  ft,  1100  gpd/ft^).  These  differences  could 
be  caused  by  the  influence  of  the  main  channel  at  their 
extremities . 


b.  Rising  and  Falling  Head  Tests:  The  rising  and  falling 
head  (slug)  tests  were  conducted  mainly  in  areas  where 
no  pump  tests  had  been  performed.  Some  testing  was 
done  west  of  the  existing  pilot  plant  and  at  other 
locations . 

The  slug  test  consists  of  placing  a  calibrated  pressure 
transducer  in  a  well  to  measure  the  water  level  in  that 
well,  removing  (or  injecting)  a  volume  of  water  from 
(into)  the  well  to  change  the  water  level  in  as  nearly 
instantaneous  a  manner  as  possible  and  recording  the 
recovery  of  the  water  level  to  its  original  value  with 
the  passage  of  time.  The  continuous  record  of  water 
level  versus  time  is  then  plotted  as  the  ratio  of 
measured  head  of  water  in  the  well  to  the  initial  head 
of  water  upon  withdrawal  (or  injection)  at  time  zero 
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(called  the  "recovery  ratio"  or  H/H  )  versus  the 
logarithm  of  elapsed  time  in  seconds.  The  curvilinear 
graph  is  then  matched  to  a  previously  calculated  family 
of  theoretical  surveys  that  includes  the  variables  of 
coefficient  storage,  transmissibility ,  permeability, 
and  confining  conditions . 

Upon  successfully  matching  the  field  data  plot  to  one 
of  the  theoretical  type  curves,  the  nature  of  aquifer 
confinement  is  identified  by  the  shape  of  the  curve  and 
the  value  of  "  "  for  the  matched  curve.  Also  the  value 
of  time  is  noted  on  the  data  plot  which  coincides  with 
the  time  of  1.0  sec  on  the  theoretical  type  curve. 


If  the  groundwater  response  in  the  aquifer  during  the 
short  period  of  the  test  and  for  the  small  volumes  involved 
indicates  unconfined  conditions,  then  the  proper  type 
curve  can  be  matched  so  that  a  value  of  the  coefficient 
of  permeability  can  be  obtained  from  the  equation; 


k  «  £n 


m  afn! 
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where 


3  2 

k  =  coefficient  of  permeability  Ch  /T/L  ) 

R  =  radius  of  influence  of  the  test  (L) 
e  * 

r  =  radius  of  well  (equal  to  radius  of  screen  in  all  of  the 
^  RMA  tests  (L) 

L  =  screen  length  (L) 

tj^  =  time  value  on  data  plot  coinciding  with  t  =  1.0  sec  (T) 


For  fully  penetrating  wells; 
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where 


H  =  height  of  stable  water  level  above  bottom  of  screen  (L) 

C  =  value  obtained  from  plotted  results  of  electrical  analog 
tests  for  a  specific  value  of  Wr^  (dimensionless) 
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For  partially  penetrating  wells; 


where 


A  +  B  £n 


(6) 


D  =  height  of  stable  water  level  above  bottom  of  aquifer  (L) 

A  =  value  obtained  from  plotted  results  of  electrical  analog 

tests  for  a  specific  value  of  L/r  (dimensionless) 

w 

B  =  value  obtained  from  plotted  results  of  electrical  analog 

tests  for  a  specific  value  of  L/r  (dimensionless)  1 

w 

If  the  groundwater  response  in  the  aquifer  during  the  short  period  of 
the  tests  and  for  the  small  volumes  involved  indicates  continued  storage 
conditions,  then  the  coefficient  of  storage  can  be  calculated  from; 


(7) 


whei;e 

5  5=  coefficient  of  storage  of  the  aquifer  (dimensionless) 

r  =  radius  of  casing  in  interval  of  water  level  fluctuation  (L) 
c 

r  =  radius  of  screen  (L) 

S  A 

o  =  value  obtained  from  type  curve  (dimensionless) 
Transmissibility  is  calculated  from; 


r  ”2 
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where  r  and  t  were  previously  defined  and  coefficient  of  permeability 
w 

is  calculated  from; 


k  =1 

^  L 


(9) 


where  T  and  L  were  previously  defined. 
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The  wells  used  for  slug  tests  were  8  inch  in  diameter 
and  were  backfilled  by  pea  gravel  subsequent  to  place¬ 
ment  of  the  piezometers  and  prior  to  sealing  with, 
cuttings.  The  piezometer  risers  and  screens  were 
1.0  inch  inside  radius.  Therefore,  in  the  analyses 
with  Equations  4,  5,  6,  7,  and  8  it  was  assumed  that 

r  =  r  =  r  =  1.0  in.  =  2.54  cm. 
c  s  w 

The  screen  lengths  were  used  for  the  values  of  L  in 
equations  4,  5,  6,  and  9  were  determined  in  the 
following  way.  A  f-ft  screen  section  was  measured  and 
found  to  consist  of  85  percent  of  its  total  length 
comprising  the  slotted  portion  and  the  remainder  of 
the  total  length  being  solid  end  sections  and  couplings 
Therefore,  the  screen  lengths  of  each  piezometer  were 
multiplied  by  0.85  to  obtain  the  value  used  for  L  in 
that  particular  calculation.  In  a  few  instances,  the 
screen  sections  extended  into  the  lower  Denver  clay 
shale  or  aquiclude.  In  those  instances  the  elevation 
differences  between  the  tops  of  the  screens  and  the 
aquiclude  were  used  as  the  nominal  screen  lengths  to 
which  the  0.85  adjustment  was  applied.  All  the  above 
assumptions  and  considerations  have  been  applied  to  all 
previous  slug  test  analyses  from  RMA  piezometers  as 
described.  This  procedure  makes  all  the  results  of  KMA 
slug  test  analyses  internally  and  directly  comparable, 
at  least  as  far  as  the  design,  construction,  and  final 
configuration  of  the  piezometers  are  comparable. 

The  bailer  used  at  RMA  for  the  slug  tests  had  a  nominal 
2.0-liter  capacity.  This  volume  of  extracted  water 
caused  an  initial  water  level  change  of  approximately 
3.2  ft  in  the  casing.  Variations  in  the  degree  of 
filling  the  bailer  coupled  with  water  level  recovery 
during  the  1-2  sec  allowed  for  surging  and  dribble- 
back  resulted  in  the  initial  water  level  being  up  to 
0.5  ft  less  than  the  masimum  possible.  The  pressure 
transducer,  together  with  the  resolution  of  the  con¬ 
tinuous  oscillographic  recorder,  provided  a  measured 
precision  of  ^0.01  ft  head  of  water.  Depending  upon 
the  time  scale  used  for  a  particular  test  the  precision 
in  time  measurements  was  either  4D. 01  or  HK).10  sec  with 
the  latter  most  commonly  used.  The  transducer  and 
recorder  were  used  for  water  level  measurements  from 
the  initiation  of  a  test  to  either  its  completion 
(judged  to  be  95  percent  recovery)  or  3000  sec  elapsed 
time,  whichever  came  first.  If  the  test  had  not  reached 
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completion  in  3000  sec,  then  the  M-scope  water  level 
detector  was  used  at  periodic  intervals  thereafter  to 
completion.  Each  initial  M-scope  reading  was  made  while 
the  transducer  data  was  still  being  recorded  to  provide 
a  consistent  date  initial  base.  It  was  found  that  the 
M-scope  data  was  reliable  to  about  +^.05  feet  of  water 
level  and  the  time  reliability  good  to  about  +1  min. 

c.  Permeability  values  of  the  north  boundary  aquifer  system 
utilized  in  the  mass  flux  computation  are  as  follows: 

_ UNIT  _ PEBMEABILITY  (FT/DAY) 

Impermeable  areas 
Silty  sands  (SM) 

Slightly  silty  sands  (SMSP) 

Sands  (SP) 

Sands  and  gravel  (SPGP) 

These  permeabilities  are  based  on  evaluation  of  the  results 
of  five  pump  tests  and  approximately  60  slug  tests.  Permea¬ 
bility  values  for  each  interval  are  summarized  in  Table  3. 

Contaminant  Concentration  Assessment 

Fluoride  in  groundwater  at  RMA  has  been  limited  to  leakage  of  fluoride 
wastes  from  disposal  basins  and  to  desorption  of  fluoride  from  fluoride 
enriched  natural  soils.  Since  the  government’s  contribution  to  an  off-post 
migration  problem  has  been  unknown,  wells  were  placed  along  the  north  boundary 
(example  of  special  distribution  is  presented  in  Figure  13) .  Routine  water 
quality  samples  from  these  wells  have  established  Basin  F  as  a  major  source 
of  the  groundwater  pollution.  Concentration  data  averaged  over  the  years 
1977  to  1979  was  used  as  input  to  the  mass  flux  computation.  A  select 
presentation  of  this  fluoride  concentration  data  along  the  north  boundary 
(Figure  14)  and  an  isoconcentration  map  for  the  entire  north  boundary 
vicinity  (Figure  15;)  are  discussed  below. 

’1 

a.  Two  areas  containing  fluoride  concentrations  in  excess  of 
5  mg/X-  are  found.  One  area  is  located  in  the  northeast 
corner  of  section  23.  These  two  areas  were  probably  connected 
at  some  time  in  the  past  and  have  been  separated  by  ground- 
water  movement.  The  concentration  found  in  the  center  of 
section  23  is  in  an  area  of  slow  flow  and  the  contaminants 
have  remained  somewhat  stationary.  The  fluoride  probably 
migrated  into  this  area  when  groundwater  elevations  were  higher. 
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b.  There  is  a  general  distribution  of  fluoride  above  2  mg/£ 
in  wells  located  along  the  northern  perimeter  of  Basin  F, 
through  most  of  section  23  and  the  western  part  of  section 
24,  to  the  northern  boundary.  The  highest  concentrations 
of  fluoride  appear  to  be  -crossing  the  transect  line  to 

the  west  of  the  existing  interim  treatment  system.  No 
well  defined  plume  was  found  exiting  the  northeast  corner 
of  Basin  F  which  is  surprising  since  Basin  F  liquid  con¬ 
tains  fluoride  in  concentrations  in  excess  of  100  mg/ 5-. 
Several  wells  in  section  24  east  of  the  main  area  of 
contaminant  distribution  were  found  to  contain  fluoride 
concentrations  above  background  levels.  Several  peri¬ 
meter  wells  on  the  southwest  corner  of  Basin  F  were 
found  to  contain  fluoride  in  excess  of  4  mg/ii  indicating 
possible  migration  in  that  direction. 

c.  Although  Basin  F  still  appears  to  be  the  source  of  fluordie 
contamination  in  study  area,  the  amount  of  fluoride  leaching 
to  the  groundwater  appears  to  be  decreasing.  The  allowable 
flux  as  indicated  on  the  flux  diagrams  will  be  shown  to 

be  very  close  to  the  actual  flux. 
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MASS  FLUX  MODELING  PROCEDURE 


roundwater  flow  through  each  aquifer  interval  is  computed  directly  from 
geohydrologic  data  provided  to  the  software  as  input.  Contaminant  concen¬ 
tration  data  is  selected  from  the  USATHAMA  Data  Management  System  according 
to  an  input  selection  routine.  If  the  model  determines  that  more  than  one 
well  exists  within  an  inteval,  a  weighted  concentration  average  is  used: 


C  = 
m 


C  N  +  C,  N, 
a  a  b  b 

N  +  N^ 
a  b 


where: 

C  =  mean  concentration  within  interval 
m 

C  ,  C,  =  concentration  at  wells  a  and  b 
a  b 

N  ,  N,  =  number  of  water  quality  values  at  well  a  and  b 
a  b 

Appropriate  multiplication  of  each  interval's  flow  calculation  with  contam¬ 
inant  concentration  data  yields  a  computed  mass  flux.  A  mass  flux  curve 


is  generated  when  all  the  intervals  are  considered.  If  the  user  specifies 
a  constant  concentration  value  equal  to  the  water  quality  standard  as  input, 
a  comparison  curve  will  be  provided.  This  comparison  flux  can  be  used  to: 

a.  Determine  locations  along  the  transect  where  the  largest 
contaminant  migration  is  likely. 

b.  Provide  an  insignt  into  the  relative  flow  of  groundwater 
across  the  transect. 


Integration  of  the  computed  mass  flux  curve  will  yield  an  estimate  of  total 
contaminant  release  through  the  section  of  aquifer  under  consideration. 
Division  of  this  estimate  by  the  groundwater  flow  rate  will  result  in  a  con¬ 
taminant  loading  expectation  expressed  in  concentration  terms. 

4.4  COMPUTATION  METHODOLOGY 

Geohydrologic  definition  cannot  be  firmly  established.  Various  individuals 
may  look  at  the  same  data  base  and  make  divergent  interpretations.  Because 
of  this,  a  validation  procedure  for  the  mass  flux  model  was  used  to  ensure 
computation  methodology  accuracy . 


An  initial  mass  flux  computation  was  conducted  using  idealized  hydrodynamic 
parameters  along  the  north  boundary.  Performance  of  the  1500  foot  pilot 
containment/treatment  system,  which  exists  within  intervals  15—21,  was 
used  for  validation  purposes.  Actual  influent  flow  rates  and  concentration 
loadings  for  the  pilot  system  were  compared  to  predicted  values  from  the  mass 
flux  computation.  Only  slight  modifications  to  the  input  parameters  were 
needed  to  achieve  an  agreement  within  a  20  percent  error.  Resultant  values 
for  the  validation  run  are  shown  below: 


EXPECTED 

FLUX 


ACTUAL 

PERFORMANCE 


Pilot  System  Influent  Flow  (gpra)  110  30-60 

Pilot  System  Influent  Fluoride  3.6  3.1-4 

Concentration  (mg/ £) 
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(1)  NOTE;  D'Appolonia  0-979)  examined  performance  of  the  pilot  contain¬ 
ment  system  over  the  first  year  of  operation.  It  was  concluded 
that  the  pilot  system  was  effectively  removing  and  treating 
approximately  30  percent  of  the  groundwater  flow  in  the  1500 
foot  aquifer  interval  intercepted  by  the  system.  An  improper  | 

match  between  aquifer  characteristics  and  dewatering  well  com-  | 

ponents  has  resulted  in  greater  groundwater  flows  than  the  \ 

pilot  dewatering  wells  could  handle- 

Once  hydrodynamic  input  parameters  were  validated,  the  mass  flux  computation 
was  utilized  to  assess  contaminant  migration  potential  and  to  perform 
various  contaminant  loading  calculations.!  Th^  following  section  expands 
on  these  assessments. 


4.5  TOTAL  QUANTITIES  OF  CONTAMINANT  TO  BE  TREATED 

Total  flow  across  the  north  boundary  has  been  estimated  from  geohydrologic 
data  presented  in  Sections  2  and  4  at  882,200  gallons  per  day  or  612  gallons 
per  minute.  This  calculation  agrees  closely  with  Zebell's  (.1979)  estimate 

*884,100  gallons  per  day.  Integration  of  the  resulting  mass  flux  computation 
Igure  16)  reveals  that  approximately  7.8  kg  of  fluoride  is  contained  in 
the  north  boundary  alluvial  flow  (Table  4) .  This  corresponds  to  an  average 
concentration  in  all  intervals  of  2.3  mg/Jl.  High  flow  rates  of  relatively 
clean  groundwater  in  the  First  Creek  vicinity  combined  with  low  flow  rates 
of  contaminated  groundwater  from  the  Basin  F  area  yield  an  expected  com¬ 
posite  stream  to  the  expanded  north  boundary  system  just  below  the  State  of 
Colorado  standard  of  2.4  mg/Jl.  Because  of  the  extensive  water  quality  and 
geohydrologic  data  available  for  this  assessment,  the  amount  of  fluoride 
to  be  processed  should  be  fairly  accurate. 


Changes  in  the  expected  contaminant  loading  to  the  expanded  treatment 
facility  over  several  years  of  operation  should  be  negligible.  Fluoride 
concentration  patterns  at  the  north  boundary  have  changed  very  little  over 
the  last  few  years  of  monitoring.  Current  contaminant  contour  lines  lie 
parallel  to  the  groundwater  flow  direction  which  will  result  in  very  little 
future  water  quality  variation.  In  fact,  upon  initiation  of  source  control 
jjj^sures ,  contaminant  loadings  at  the  north  boundary  expanded  system  will 
significantly  lessen. 
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TABLE  3 

FLUORIDE  MASS  FLUX  INPUT  PARAlffiTERS 


Transect  Length  10,000  Feet 

Interval  Width,  250  Feet 


Water  Quality  Sampling 

Interval 

1  Jan  77 

to  4  Oct  79 

Water  Quality  Standard 

2.4  mg/£ 

INTERVAL 

PERMEABILITY 

GRADIENT 

THICKNESS 

LENGTH 

NO. 

(FT/DAY) 

CFT/FT) 

(FT) 

■  (FT) 

1 

10.0 

0 

0 

0 

2 

10.0 

0 

0 

0 

3 

10.0 

0 

0 

0 

4 

10.0 

0 

0 

0 

5 

10.0 

0 

0 

0 

6 

10.0 

0 

0 

0 

7 

0.4 

0 

0 

0 

8 

0.4 

0 

0 

0 

9 

0.4 

0 

0 

0 

10 

50 

0 

0 

0 

11 

0.4 

.0030 

1.0 

100 

12 

50 

.0040 

3.0 

250 

13 

50 

.0050 

3.0 

250 

14 

100 

.0060 

7.0 

250 

15 

150 

.0062 

10.0 

250 

16 

150 

.0064 

12.0 

250 

17 

200 

.0066 

11.0 

250 

18 

250 

.0068 

9.0 

250 

19 

300 

.0070 

8.0 

250 

20 

350 

.0072 

7.0 

250 

21 

350 

.0074 

7.0 

250 

22 

300 

.0076 

8.0 

250 

23 

225 

.0078 

10.0 

250 

24 

350 

.0080 

12.0 

250 

25 

225 

.0078 

16.0 

250 

26 

180 

.0076 

17.0 

250 

27 

250 

.0074 

12.0 

250 

28 

250 

.0072 

15.0 

250 

29 

200 

.0070 

18.0 

250 

30 

225 

.0068 

21.0 

250 

31 

350 

.0066 

22.0 

250 

32 

350 

.0064 

20.0 

250 

33 

325 

.0062 

15.0 

250 

34 

250 

.0060 

6.0 

250 

S5 

200 

.0058 

2..0 

50 

36 

0.4 

0 

0 

0 

37 

0.4 

0 

■  0 

0 

38 

0.4 

0 

0 

0 

39 

0.4 

0 

0 

0 

40 

0.4 

0 

0 

0 
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TABLE  4 


TOTAL  QUANTITIES  OF  FLUORIDE  TO  BE  PROCESSED 


INTERVAL 

flow' RATE  1 

FLUORIDE  MASS  FLUX 

NO. 

(<3PD) 

(GRAMS /DAY) 

1 

0 

0 

2 

0 

0 

3 

0 

0 

4 

0 

0 

5 

0 

0 

6 

0 

0 

,7 

0 

0 

8 

0 

0 

9 

0 

0 

10 

0 

0 

11 

1 

0 

12 

1120 

tc - 

13 

1400 

50 

14 

7850 

90 

15 

17400 

280 

16 

21550 

340 

27140 

350 

fl8 

28500 

390 

19 

31400 

440 

20 

33000 

340 

21 

33910 

420 

22 

34110 

370 

23 

32820 

300 

24 

62830 

650 

25 

52500 

380 

26 

43380 

410 

27 

42630 

440 

28 

50490 

450 

29 

47120 

420 

30 

60100 

400 

31 

95030 

490 

32 

83650 

370 

33 

56550 

275 

34 

16850 

125 

35 

869 

10 

36 

0 

0 

37 

0 

0 

38 

0 

0 

39 

0 

0 

40 

0 

0 

TOTAL 


882200 


7790 


INTERVAL 


tonm  MMtRS 


RMA  NORTH  BOUNDARY 


•  w 
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FIGURE  13  RMA  NORTH  BOUNDARY  MONITORING  WELLS 


RMA  NORTH  BOUNDARY 
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FIGURE  14  MEAN  FLUORIDE  CONCENTRATION  (mg/I)  AT  NORTH  BOUNDARY  RMA  (1977-1978) 


MASS  FLUX  OF  FLUORIDE  ACROSS  NORTH  BOUNDARY  OF  RMA 


(nvAUiXNi/Ava/siAio)  xmj  ssvw  aaiaomd 
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5.0  CONCLUSIONS 

A.  Groimdwater  flow  at  BMA  is,  predominantly  from  south  to  north.  Locally 
in  the  north  boundary  vicinity  two  separate  subsurface  water  units  make  up 
the  alluvial  flow  crossing  the  Arsenal’s  northern  boundary.  One  of  these 
units  moves  beneath  Basin  F  in  a  northeasterly  direction.  Contaminants 
leached  from  surface  waste  basins  move  within  this  flow.  The  other  groundwater 
unit  is  relatively  free  of  contaminants  as  it  proceeds  parallel  to  First  Creek. 
Therefore,  water  quality  within  the  alluvial  aquifer  at  the  boundary  is  greatly 
dependent  on  the  pathways  taken  by  the  individual  groundwater  units. 


B.  Total  north  boundary  alluvial  flow  control  is  envisioned  to  meet  applicable 
State  of  Colorado  water  quality  guidelines.  Interception  of  the  entire 
alluvial  aquifer  will  result  in  compositing  both  the  contaminated  and  relatively 
noncontaminated  groundwater  flows  units  described  above.  Independent  assess- 
^tets  have  been  completed  for  design  purpose  to  predict  expected  flow  and 
^Ritaminant  loading  within  the  expanded  north  boundary  control  scheme. 


C.  Flow  estimates  of  the  alluvial  aquifer  at  the  north  boundary  are  between 
450  and  650  gpm.  Variation  of  these  estimates  ia  due  primarily  to  the  choice 
of  hydrodynamic  parameters  for  the  aquifer.  Permeability  estimates  for  the 
most  permeable  aquifer  material  range  from  400  to  600  feet  per  day .  Equal 
variation  is  noted  in  saturated  thickness.  These  differences  are  within 
reason,  however.  Geohydrologic  definition  is  not  an  exact  science  and  is 
assumed  adequate  if  an  80  percent  accuracy  is  achieved. 


D.  Total  mass  of  fluoride  contained  in  the  alluvial  aquifer  as  it  passes 
off  RMA  is  estimated  at  6.6  to  7.8  kg  per  day.  This  equates  to  an  average 
fluoride  water  quality  of  2.3  to  2.7  mg/a.  Because  of  the  geohydrologic 
variation  previously  noted,  refinement  of  these  concentration  expectations 
will  not  be  possible  until  actual  expanded  system  operation  is  accomplished. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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APPENDIX  A 

FINITE  ELEMENT  PROGRAM  DESCRIPTION 
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NOMENCLATURE 


SYMBOL 


MAXBW 


DESCRIPTION  DIMENSION 

Coefficients  in  Equation  (1.5.14)) 

Thickness  of  aquifer  ^ 

Thickness  of  semipervious  layer  ^ 

Domain  of  interest 

Superscript  representing  an  element 

Exponential  integral 

Subscript 

L/t 

Accretion 

Subscript 

Subscript 

Principal  hydraulic  conductivity  along  x  axis  L/t 

Principal  hydraulic  conductivity  along  y  axis  L/t 

Hydraulic  conductivity  of  the  semipervious  layer  L/t 
Directional  cosine  with  respect  to  x 

Directional  cosine  with  respect  to  y 

L 

Length 
An  integer 

Number  of  nodes  in  one  element 

Order  of  approximation 
•  Upper  bandwidth  of  matrix  plus  one 

Subscript 
Shape  function 

Number  of  elements  in  whole  grid  system 
Number  of  nodes  in  whole  grid  system 
Shape  function  for  a  specific  point 

2 

M/Lt^ 

Pressure  -a  o 

L  /L  t 

Discharge  per  unit  area  ^  2 

Dischatge  par  unit  area  at  a  specific  point  I-  /i  t 


Known  flux  along  boundary 


1.0  PROGRAM  INFORI'IATION 


1  1  npTrrw  A\rn  purpose  OF  PROGR^ 

computer  progrem  by  the  code  name  FICA  (Floo  in  Confined  Aqurfer)  van 
developed  at  Michigan  State  Dnlversity.  Department  of  Cxvrl  Engrneerrng 
by  Sirous  Baji-Djafari  and  David  C.  Wiggert.  This  program  is  a  modi,  le 
and  improved  version  of  the  already  documented  program  by  David  C. 
Kiggert.  This  version  uses  isoparametric  elements  (quadrilateral,  tri 
angular  or  mined  elementsl  while  the  previous  one  only  employed  linear 
triangular  elements.  Some  of  the  subroutines  of  the  program  are  provi  e 
by  Dr.  L.  J.  Segerlind  of  the  Department  of  Agricultural  Engineering. 

Michigan  State  University. 

The  purpose  of  the  program  is  to  simulate  the  performance  of  an 
aquifer  on  a  regional  basis  with  a  two-dimensional  model.  Finite  ele¬ 
ment  method  is  used  to  discretize  the  governing  partial  differential 
equations.  By  providing  hydrodynamic  parameters  and  stresses  (such  as 
transmissibility.  storage  coefficient,  pumping  rate,  etc.),  the  program 
will  find  piezometric  heads  at  the  nodes  and  consequently  the  velocity 
vectors  either  at  the  nodes  or  within  the  elements.  In  addition,  a 
steady-state  solution  can  be  computed  either  as  an  initial  or  final 

condition. 


1,2  A-REAS  OF  APPLICATION 

The  program  FI*CA  can  be  used  to  - simulate  two-dimensional  groundwater 
flow  in  anisotropic  and  nonhomo geneous  aquifers  under  confined  or  un¬ 
confined  aquifers.  Leaky  artesian  aquifers  can  be  incorporated  in  the 
program.  The  other  features  of  the  program  include  time  variable 
p-ompage  from  well,  natural  or  artificial  recharge  and  line  source 
recharge, 

Tne  Galerkin  finite  element  formulation  is  employed  to  discr 
space  and  time  derivatives  of  governing  equations.  Any  isop 
element  can  be  used  for  grid  system  development. 
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1,3  PROBLEM  DESCRIPTION  AND  USER  ORIENTATION 

mathematical  equations  describing  the  flow  through  a  confined  or  un- 
T^fined  aquifer  in  two-dimensional  horizontal  planes  are  described  in 
Section  l.A.  Finite  element  formation  of  these  equations  are  .developed 
in  Section  1.5.  The  computer  model  solves  the  flow  equation  in  the 

following  procedure. 


The  domain  of  investigation  is  divided  into  a 
group  of  isoparametric  elements.  These  elements 
can  be  triangular,  quadrilateral  or 

cents  of 'any  kind  (i.e.,  linear,  quadratic  or  cubic) 
By  means  of  these  elements  and  employing  the 
Galerkin-based  finite  element  method,  the  mathe¬ 
matical  equations  are  transformed  to  a  system 
of  first-order  partial  differential  equations. 

The  variation  <})  within  the  element  depends 
on  the  kind  of  element.  Other  properties  such 
I  as  transmissibility,  storage  coef f icient ,  and 
'  recharge  are  assumed  constant  within  an  element. 

To  discretize  the  recurrence  formula,  different 
techniques  are  discussed  (see  Equation  1.5.1  ). 

The  user  has  the  option  to  choose  any  of  the 
presented  methods. 


Some  of  the  features  of  the  program  are  as  follow: 

•  At  any  designated  node,  time-vafiable  pumping 
rates  can  be  imposed. 

•  The  system  can  be  divided  into  regions  and  each 
region  can  have  different  known  parameters  such 
as  transmissibility,  recharge  and  storage 
coefficient. 

•  Initial  piezometric  head  can  be  specified  at  all 
nodes  or  can  be  computed  by  the  program. 


At  any  designated  node,  the  value  of  the  piezo¬ 
metric  head  can  be  specified  (Dirichlet  boundary 

condition) . 

Along  any  designated  element,  line  source  can  be 
incorporated. 

Piezometric  heads  can  be  computed  either  for  a 
steady-state  condition  or  for  a  transient  state. 
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•  Velocity  vectors  can  be  computed  at  any  specified 
time  either  at  the  nodes  or  in  the  elements. 

Basic  assumptions  used  for  developing  this  program  are  presented  in 
Section  1.4.1.  Among  those,  it  is  assumed  that  the  variation  of  thick¬ 
ness  of  the  aquifer  is  much  smaller  than  the  thickness  itself.  Further¬ 
more,  the  vertical  flow  components  are  of  minor  importance  and  flow 
components  are  essentially  two  dimensional  and  horizontal. 

When  predicting  drawdown  close  to  the  wells,  true  drawdown  cannot  be 
accurately  computed  at  the  well  node.  In  order  t(5.  obtain  greater  accuracy , 
the  smaller  nodes  are  used  in  the  vicinity  of  the  well  point. 

In  the  case  of  phreatic  flow,  it  is  important  to  note  that  the  related 
equations  are  approximate  only  and  their  use  may  result  in  significant 
error  if  computed  drawdown  becomes  large  relative  to  the  initial  saturated 
thickness.  The  present  program  can  be  modified  to  permit  the  calculation 
of  the  apparent  transmissibility  defined  by  Equation  1.4.8.  In  this 
case,  the  average  piezometric  head  is  calculated  at  every  time  step  and 
multipicked  by  hydraulic  conductivity.  If  more  accuracy  is  desired,  it 
is  possible  to  iterate  within  one  time  step,  each  iteration  producing  an 
update  value  for  4),  which  is  used  to  reevaluate  the  apparent  trans¬ 
missibility.  This  procedure  requires  regenerating  and  decomposing 
global  matrices  when  new  transmissibilities  are  computed,  and  hence, 
increasing  processing  time. 

Section  2.0  of  this  manual  deals  with  usage  information.  To  simplify 
input  data  preparation,  two  tables  are  presented.  Table  II-l  shows  the 
name,  location  and  order  of  variables  along  with  the  number  of  data 
cards  and  format  numbers.  The  required  format  for  each  order  is  shown 
in  Table  11-2.  A  sample  problem  is  given  in  more  detail  to  orient  the 
user  with  the  steps  which  are  required  to  use  this  program. 

Section  3.0  contains  the  listing  of  the  programs- 
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l.A  BASIC  GOVERNING  EQUATIONS  OF  GROUNDVATER  FLOW 

regional  problems,  two-dimensional  horizontal  flow  is  considered. 
The  governing  equations  are  well  established  (e.g.,  see  Bear,  1972,  and 
Pinder  and  Frind,  1972). 


1.4.1  Basic  Assumptions 

The  following  assumptions  are  valid  for  regional  groundwater  flow: 

(a)  The  flow  is  essentially  horizontal  in  a 
two-dimensional  plane  (Figure  I-l).  This 
assumption  is  valid  when  the  variation  of 
thickness  of  the  aquifer  is  much ^smaller  than 
the  thickness  itself.  This  approximation 
fails  in  regions  where  the  flow  has  a  vertical 
component. 

(b)  The  fluid  is  homogeneous  and  slightly 
compressible. 

(c)  The  aquifer  is  elastic  and  generally  non- 
homogeneous  and  anisotropic.  The  consolidating 
medium  deforms  during  flow  due  to  changes  in 
effective  stress  with  only  vertical  compress¬ 
ibility  being  considered. 

(d)  For  the  two-dimensional  horizontal  flow 
assumption,  an  average  piezometric  head  is  used 
where  the  average  is  talcen  along  a  vertical 
line  extending  from  the  bottom  to  the  top  of 
the  aquifer,  i.e.. 


4><x,y,z,t)dz 


(1.4.1) 


where  b  is  the  thickness  of  the  aquifer. 


The  piezometric  head  is  defined 

({)  =  ^  +  Z 
Y 


where  p  is  pressure,  y  unit  weight  of  fluid,  and  Z  elevation  from  a  datum. 
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1.4.2  T\s’o-Dinienslonal  Horizontal  Flow  -  Confined  Aqu^  _ 

The  combined  equation  of  motion  and  continuity  for  flow  in  a  two 

dimensional  horizontal  plane  can  be  written 


(1.4.2) 


„here  K  ,  K  .re  principal  hyaraulic  conductivities  along  x  and  y 
directirn.  /L  thiCness  of  confined  aquifer,  P  is  the  s^engt  o 
einh  (or  source,.  1  is  the  vertical  recharge 

aquifer.  is  the  elastic  of 

derivatives  with  respect  to  x,  y»  »  ^  -.-i  j  t.v 

hydraulic  conductivity  and  the  thickness  of  aquifer  is  ca 

transmissibility  of  the  aquifer.  Thus 


T  =  bK  and  T  -  bK 

XX  XX  yy  >> 


(1.4.3) 


and  storage  coefficient  is  defined  by 


S  =  S  b 
s 


(1.4.4) 


For  a 


confined  aquifer.  Equation  (1.4.2)  becomes 

ft 

3  /■.p  +  — (T  1^)  -  P  +  I  = 

I^^^xx  3x^  3y^  yy  9y 


(1.4.5) 


1.4,3  Initial  and  Boundary  Conditionj_ 


Joundary  Conditions 

In  order  to  solve  a  partial  differential  equation  describing  a  p  ys 

choose  certain  additional  conditions  imposed 
phenomenon,  it  is  necessary  to  choose  certai 

hy  the  physical  situation  at  the  boundaries  (S)  for  the  d,^  n  0)  und 
uLsideration.  In  general  the  equation  for  the  boundary  condrtron 

written 
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(T  +  T  +  6_<})  +  B  =  0 

1  XX  dx  i  yy  oy  2  2  3 


(1.4.6) 


^1  ’  ^2  tl'.e  directions!  cosines,  snd  B. »  B.,,  snd  B^  sre  given 
functions  of  position  and  possibly  time.  For  flow  through  an  aquifer, 
thrc'e  different  boundary  conditions  are  applicable: 

(a)  Dirichlet  or  prescribed  potential:  In  this 

case  the  potential  is  specified  for  all  points 
along  the  boundary 
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(b)  Neumann  or  prescribed  flux:  Along  a  boundary 
•^f  this  type,  the  flux  normal  to  the  boundary 
surface  is  prescribed  for  all  points  of  the 
boundary  as  a  function  of  position  and  time 


T  +  T  = 

XX  9x  1  yy  3y  2 


-  on  S;  B,  0 


A  special  case  of  the  Neumann  condition  is  the 
impervious  boundary  where  the  flux  vanishes 
everywhere  on  the  boundary,  i.e.. 


(c)  Cauchy  boundary:  This  problem  occurs  when  the 

potential  and  its  normal  derivative  are  prescribed 
on  the  boundary  in  the  combined  form,  and  the 
entire  Equation  (1.4.6)  is  used. 

forms  of  Equation  (1.4.6)  for  three  types  of  boundary  conditions 
are  summarized  in  Table  I-l.  In  general,  for  a  flow  problem  one  will  have 
mixed  boundary  conditions  in  which  the  Dirichlet  condition  will  apply  over 
a  part  of  the  boundary  and  the  Neumann  condition  will  be  specified  for  the 
remaining  portion  (Bear,  1972). 
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either  the  piatot,etrlc  heeds  ere  keo™  xn  the  ee  tre 

do^it  (D)  or  the  hPdroloslc  stresses  (seeh  es  poppies 
specified  and  boundary  conditions  are  known.  For  t  e  sec 
system,  has  reached  the  steady  state,  so  the  solution  of  t  e  equa 


li)  +  ^(T  |i)  -  P  +  I  =  0 

3x^  XX  3x^  3y  yy  3y 


(1.4.7) 


will  yield  piezometric  heads  for  the  initial  time. 

1  A  4  Tt.^-T^^’^^n.lonal  Horizontal  Flow  -  Uncojifitied_AaHiiil 
T^e'aquI^ITirown  in  Figure  X-2  is  for  saturated  flow  bounded  above 

3  phreatic  surface.  For  regional  analysis  it  is  possible  to  descri  e 
1  flow  With  a  relation  analogous  to  Equation  (1.4.2)  by  making  use  of 
the  Dupuit  approximation.  In  this  case  the  transmissibility  becomes 


T  =  4>K 

XX  XX 


T  =  4>K 

yy  yy 


(1.4.8) 


vhere  it  is  assumed  the  impervious  boundary  is  the  datum.  In  Equation 
(14  8)  *  is  the  pietometric  head  at  any  specified  location,  i^en 

ap;i;ii;  !  .nss  balance  to  a  control  volume  in  Eisure  1-2.  in  addition 
to  compressibility  of  the  fluid  and  porous  media,  the  variation  of  ava  l- 
able  storage  due  to  vertical  movement  of  the  phreatic  surface  should  be 
considered.  The  concept  of  drainable  mater  or  equivalently  specific  yi  . 
S  .  can  be  used  to  describe  this  phenomenon.  Since  in  most  cases  »  S. 
in  Equation  (1.4.2)  can  be  substituted  for  S. 


.  Equation  (1.4.2)  takes  the  form  of 


S/v-  +  — (K  <i^)  “  ^  ^  ft  (1.4.9) 
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^fundamental  difference  occurs,  however,  due  to  the  dependence  of  trans¬ 
missivity  upon  4*  3S  shown  in  Equation  (1.4.8).  Thus  the  unconfined  flow 
equation  is  nonlinear,  and  numerically  is  treated  in  a  quasi-linear 
fashion.  In  unconfined  conditions,  infiltration,  I,  directly  enters  the 
saturated  zone.  Infiltration  can  be  either  natural  or  artificial  input 

to  the  system. 


1.4.5  Infiltration 

Unconfined  aquifer  recharge,  I,  directly  infiltrates  to  the  aquifer.  In 
nonleaky-confined  aquifers,  the  recharge  value  is  zero.  In  semiconfined 
(leaky)  aquifer  (see  Figure  1-3) ,  the  vertical  leakage  is  calculated 
(Bear,  1972)  by  the  Equation 


4,0-4) 


(1.4.10) 


where ; 


=  potential  head  in  the  aquifer  above  the 
semipervious  layer 

=  piezometric  head  in  confined  aquifer 

=  the  resistance  of  the  semipervious  layer 
and  is  equal  to  b’/k’ 

=  thickness  of  semipervious  layer 

=  hydraulic  conductivity  of  the  semipervious 
layer 


1.4.6  Velocity  Vectors 

The  equation  of  apparent  velocity  (flow  per  unit  area)  along  x-  and  y- 
directions  can  be  written  as  follows: 


-K 

XX  3x 


(1.4.11a) 


-K 

yy  3y 


(1.4.11b) 


where 


V^  =  velocity  vector  along  x-direction 
Vy  =  velocity  vector  along  y-direction. 


The  other  terms  are  defined  previously. 

The  magnitude  and  direction  of  the  velocity  vectors  (Figure  1-4)  at  any 
given  location  and  time  can  be  obtained  by 


V  =  .  /v^  + 

V  ^  y 


a  =  arc tan  /V^  /  V  \ 


(1.4.12) 


(1.4.13) 


where 


V  =  magnitude  of  the  velocity  always  positive 

ct  =  direction  of  the  velocity  measured  from  positive  x. 


1.5  i-lETHOD  OF  COMPUTATION 

The  finite  elenent  method  is  a  numerical  technique  vhich  is. used  to 
approximate  a  continuous  partial  differential  equation  in  a  given  domain 
B  Lh  specified  boundary  conditions  along  boundaries  S.  The  hey  features 
of  the  finite  element  concept  are  (Sortie  and  de  Vries,  1973): 

1.  The  domain  is  divided  into  subdomains  or 
finite  elements,  usually  of  the  same  order. 

2.  The  trial  solution  is  prescribed  (function-, 
ally)  over  the  domain  in  a  piecewise  fashion, 
element  by  element. 

A  detailed  formulation  of  the  finite  element  method  is  given  by  Zienkievlcr 
(1971),  Notrle  and  de  Vries  (1973).  This  technique  has  been  utilired  by 
several  investigators  (Javandel  and  Witherspoon,  1971;  Finder  and  Frind 
1972;  Neuman  and  Witherspoon,  1971;  Besai,  1972;  Cheng  and  Li,  1973;  and 
France,  1971,  1974)  to  solve  transient  flow  problems  In  a  confined  or 

unconfined  aquifer. 

in  this  section  a  brief  discussion  of  the  C.lerhin  based  finite  element 
technique  is  given  and  the  method  is  used  to  discretize  the  space  deriva¬ 
tives  of  the  flow  equation.  The  simultaneous  solution  of  velocity 
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vectors  is  also  described,  i.e.,  the  Galerkin  formulation  of  the  Darcy 
-^j^s  constructed  and  velocity  components  are  calculated  at  the  nodes. 

1.5.1  The  Galerkin  Finite  Element  Method 

In  the  finite  element  technique  the  domain  D  is  divided  into  subdomains 
D®  which  are  called  elements.  Each  element  is  designated  by  nodes.  In 
this  documentation  NELS  represents  the  number  of  elements,  M  is  the  number 
of  nodes  in  each  element,  and  NNDS  stands  for  the  total  number  of  nodes 
in  domain  D  (see  Figure  1-7) . 

Consider  a  problem  of  solving  approximately  a  set  of  differential  equations 
in  nhich  the  unknomi  function  {  has  to  be  satisfied  in  the  domain  B  mith 
the  boundary  conditions  specified  along  S.  The  governing  equation  can  be 

written 


f({!|)})  =  0 

the  trial  solution  for  this  equation  be  4> 

(j)  =  [N]  f  c})}  =  I 
n=l 


(1.5.1) 


where  [N]  =  [N(x,y)]  are  shape  functions  (prescribed  functions  of  coordinates) 
and  {4)}  ={({)(t)}.is  a  set  of  M  unknown  parameters.  In  general,  the  equation 
of  residual  (or  .error)  is  formed  in  the  following  way: 


=  f  ({({.})-  f  A  4>})  =  -f  ^  0 


(1.5.2) 


The  best  solution  will  be  one  in  which  the  residual  R  has  the  least  value  ^ 
at  all  points  in  the  domain  D®.  An  obvious  way  to  achieve  this  (Zienkiewicz. 
.1971)  is  to  make  use  of  the  fact  that  if  R  is  Identically  zero  elsewhere. 


then 


f  W  R  dD  =  0 
^  -6 
D 


(1.5.3) 
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.here  W  is  any  fnnction  of  the  coordinates.  If  the 

parameters  (  .)  is  h™S  and  hEIS  linearly  independent  fnnctro  V  a  e 
chosen,  one  can  vrite  a  suitahle  number  of  simultaneous  e,uatrons 


.  J  ^  R  dD  =  i  g  W^f((Nl[4>))dD  =i0> 

Ic  —  Ij  •••jM 


(1.5. A) 


D 


D' 


•  c  f-tn-n  Tf  the  shape  function  N,  is 
V  a  P’  -it;  called  the  weighting  function- 

r  rr:  -  --^.tins  funotion.  the  process  is  termed  the  Caler  an 
procedure  uhich  is  used  henceforth.  Ihe  element  equations  can  he  assemhied  hy 


NELS  .  in 

I  /  e  «  ®  I  ■  ^ 

e=l  * 


(1.5.5) 


to  y 


vield  the  global  relations  for  domain  D. 


tS2  Finite  ElementJForinu^ 

T;.e  residual  equation  for  flow  in  a  confined  horirontal  aquifer 
(Equation  1.A.5)  with  no  vertical  recharge  can  be  written  as 


K=S||-- 


U  lit+l-d  |i)  (1.5.6) 

3x^  3y  yy  9y  i . 

numerical  approximation  of  (}..  Substituting 


The  symbol  "  represents  the 
Equation  (1.5.6)  into  Equation  (1.5  «.  one  obtains 


li 

3t 


|i  +1-  T 


+  ^  "fe  ^xx  3x1  8y  \  yy  3y 


)J 


■M  dD  =  0  k—1 » .  • 

k 


(1.5.7) 


By  use  of  the  Green  theorem. 


the  third  term  can  be  modified 


f  li  r  M  +  #1  XudD  -  -Le 

D®  3x*  3y  yy  3y  i  ^  ■> 


^XX  3x  3x 


(1.5.8) 


-  ^yy  t  h  ”yy  h) 
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The  last  tenu  in  Equation  (1.5.8)  is  nonzero  only  for  elements  which 
obtain  the  Neumann  flux  boundary  condition 


_f  N  T  1^  £  +  T  SI  e  \  ^2  (1*5.9) 

■’  e  k  \  XX  3x  1  yy  3y  2  e  x  z 


where  Q  is  kno;*Ti  flux  along  the  boundary.  Substituting  Equation  (1.5.8) 
and  Equation  (1.5.1)  into  Equation  (1.5.7)  and  raarianging  the  tert^.  one 

obtains 


dD  + 


e  '^n  1  ^xx  3x  3x  ^yy  3y  3y 


(1.5.10) 


f  e  =  W  t  ■“=  ^  V  ^  ^2 


dS  =  0 


Sinee  (.  and  its  time  darivativea  are  independent  of  the  coordinatea, 

can  be  tahen  out  of  the  integrals.  Equation  (1.5.10)  can  be  written 


matrix  form 


(1.5.11) 


where 


3N,  3N 


3N,  3N  - 
•  k  n’ 


,  p  oi,  ua.  -  _ n, 

=  i  e  rxx  dl'  3y"  '^yy  3y  3y  J 

D 


(1.5.12a) 


k,n=l, > 


fHl®  =  =  f  S  N  N  dD 

"kn  Jjje  k  n 


(1.5.12b) 


i,j=1.2 


f  -  K  ‘  \  ‘=2  ’’  \ 


(1.5.12c) 
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It  is  assumed  that  the  storage  coefficient  is  constant  throughout  the 
element  and  that  the  element  coordinate  axes  coincide  with  the  principal 
direction  of  the  transmissivity  tensor;  the  transmissivity  can  be 
defined  either  at  the  nodes  or  at  each  element.  Evaluation  of  Equation 
(1.5.12)  for  different  types  of  elements  is  discussed  by  Zienkiewicz 
(1971)  and  is  presented  in  more  detail  by  Haji-Djafari  (1976).  Upon 
evaluation  of  Equation  (1.5.12)  for  all  elements  and  transformation  to 
a  global  coordinate  system,  they  are  assembled  by  virtue  of  Equation 
(1.5.5)  into  a  global  relationship 

[B]  {(J)}  +  [H]  {||}  =  {F}  (1.5.13) 


The  parameter  {<!>}>  matrices  [B]  and  [H],  and  force  vector  {f}  are  the 
summation  of  the  corresponding  terms  in  Equation  (1.5.12)  over  all  the 
elements  in  the  Domain  D.  The  matrices  [B]  and  [H]  are  banded  symmetric. 
Equation  (1.5.13)  is  a  set  of  first  order  linear  differential  equations 
with  unknowns  {<})}  and  can  be  solved  simultaneously  at  the  given  nodes  in 
the  space  domain. 


The  recurrence  formula  for  Equation  (1.5.13)  has  the  form  (Haji-Djafari, 
1976) 


I 

\ 


[B]  + 


!i2 

At 


/ 


{(fi(t+At)}  = 


r- 


[B]  + 


f 

/ 


{tf’(t)}  +  a^2  {F(t+At)}  +  a^^  {F(t)}  + 


(1.5.14) 


In  Equation  (1.5.14),  At  is  time  step  and  a*s  are  coefficients  which 
their  value  are  given  Section  2.2-Order  12. 

^•^•3  Finite  Element  Computation  of  Velocity  Vectors 

Mathematical  equations  of  velocity  vectors  are  presented  in  Section:  1.4.6. 
^ce  the  piezometric  heads  have  been  determined,  velocity  vectors  (or 
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flow  per  unit  area)  can  be  evaluated.  Two  techniques  are  employed  to 
compute  the  velocity  vectors  by  finite  element  method. 


The  first  technique  which  is  called  "direct  method,"  piezometric  head  is 
approximated  by  Equation  1.5.1.  The  resulting  equations  will  be: 


V  =  -K  -w- 
X  XX  dx  n 


V  =  -K  -r-  4> 

y  yy  3y  n 


n  =  1,  .  .  .  M,  number 

of  nodes  in  an  element 


(1.5.15) 


[ji  £qu3rion  (1.5.15)  s  are  piezometnc  heads  at  the  nodes,  and  are 
SN  3N 

known.  The  terms  — and  — are  first  derivatives  of  shape  functions 


and  are  evaluated  at  the  point  of  interest,  usually  at  the  center  of  element. 

In  the  second  technique  which  is  termed  "simultaneous  method,"  the  Galerkin- 
based  finite  element  formulation  of  Equation  (l.A.ll)  is  developed.  This 
procedure  yields  a  set  of  equations  which  are  solved  simultaneously  to  find 
the  velocity  vectors  at  the  nodes. 

The  detailed  procedure  is  given  by  Haj i-Dj af ari ,  1976  and  the  results  are 
summarized  below. 


The  element  equations  for  x— component  of  velocity  vector  have  the  form: 


•.v‘^1  =  -F  :• 


(1.5.16) 


where 


[H]®  = 


r  N,  N  d 

Jd" 


K,  n  =  1,  2,  .  .  .  M 


(1.5.17) 


‘  =  -''x,  be  \ 

I  -7  D 


(1.5.18) 


K,  n  =  1,  2,  .  .  .  M 
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Constructing  a  global  matrix  yields 


In  Equation  (1.5.19),  IH]  is  landed  symmetric  matrix,  and  ■  .  is  known 

column  force.  The  solution  of  Equation  (1.5.19)  yields  the  x-component 
of  velocity  at  each  node  simultaneously. 

Similarly,  the  finite  element  equations  of  the  velocity  vector  along 
y-direction  have  the  forms  similar  to  Equation  (1.5.16)  through 
Equation  (1.5.19),  except  subscript  x  is  replaced  by  y. 


1.6  ACCURACY,  LIMITATIQKS  AND  RESTRICTIONS 

As  presented  in  Chapter  2.0,  the  solutions  obtained  from  this  program  are 
compared  with  the  exact  solution  by  Theis  and  good  agreement  is  found. 
However,  this  program  is  a  numerical  model  which  is  developed  using  finite 
element  method.  As  with  any  other  numerical  technique,  there  are  errors 
which  are  associated  with  size  and  kind  of  elements;  numerical  integra¬ 
tion  method,  and  discretization  of  the  recurrence  formula;  size  of  time 
step;  etc.  Experience  shows  the  accuracy  of  the  results  will  improve  by 
reducing  the  size  of^lements  and  time  stee_.  In  general,  the  program  has 
been  verified  and  the  accuracy  of  the  results  are  within  an  acceptable  range 

Limitations  and  restrictions  of  the  program  are  described  in  appropriate 
sections,  the  major  ones  being  as  follows; 

•  Simulation  of  the  aquifer  performance  is  limited 
to  regional  situations  as  long  as  the  validity  of 
basic  assumptions  described  in  Section  l.A.l  are 
maintained. 

•  Although  any  kind  of  elements  can  be  used,  their 
improper  combination  is  restricted  as  described 
in  Section  2.2. 

•  Units  of  the  input  parameters  should  be  consistent. 
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1.7  PROGRAM  ORGANIZATION 


The  program  consists  of  main  program  PICA  and  15  internal  subroutines. 
Main  program  functions  as  an  organizator;  it  reads  job  specification 
parameters;  sets  the  dimensions  of  the  variables;  calls  for  input  sub¬ 
routine;  initializes  and  increments  time;  calls  for  subroutines  such  as 
those  identifies  matrix  bandwidth,  -constructs  global  matrices,  constructs 
force  vector,  solves  the  piezometric  heads,  prints  the  piezometric  heads, 
calculates  the  discharge  vectors;  and  finally,  the  main  program 
terminates  the  job  if  there  are  no  new  data  set  for  computation.  The 
flow  chart  for  the  main  program  is  given  in  Figure  1-5. 

The  function  of  each  subroutine  is  described  in  the  computer  program 

(see  Section  3.0  for  program  listing). 

•  : 

1.8  CASE  STUDIES 

1.8.1  Regional  Aquifer  System 

The  domain  illustrated  in  Figure  1-6  is  a  confined  aquifer  approximately 
18,000  m  in  length  and  5000  m  wide.  A  well  field  is  located  at  E,  which 
pumps  18,930  m^/d  from  the  system.  Aquifer  properties  are  K  =  50  m/d, 

S  =  0.001,  and  two  thicknesses  of  10  and  20  m,  resulting  ih  two  trans- 
missibility  zones  of  500  and  1000  m^/d,  respectively.  No  flow  boundaries 
exist  along  AjJ'anci  CD,  and  known  potential  conditions  are  assigned  along 
AD  and  BC.  The  system  is  divided  into  isoparametric  linear  elements 
yielding  Ilf  nodes  and  lOA  elements  (Figure  1-7).  The  predicted  drawdown 
at  two  locations  at  points  G  and  F  for  40  days  after  pumping  is  shown 
in  Figure  1-8. 

The  same  aquifer  is  simulated  using  the  method  of  characteristics 
(Wiggert  &  Wylie,  1976).  The  results  obtained  from  the  finite  element 
method  are  used  to  justify  the  capability  of  the  method  of  character- 
is'tics.  Good  agreement  is  found  between  the  results  of  these  two  models, 
as  well  as  an  analytical  solution,  as  shown  in  Figure  1-8. 
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1.8.2  Single  Horizontal  Drain 

The  second  case  study  involves  the  simulation  of  flow  movement  to  a  single 
horizontal  drain  18.3  m  in  length.  At  the  initial  time,  the  water  level 
is  assumed  to  be  24. A  m  above  the  barrier  boundary.  Water  is  allowed  to 
discharge  from  only  seven  points  along  the  drain,  as  shown  in  Figure  1-9. 
Known  potential  conditions  are  assigned  along  AB. 

The  system  is  divided  into  isoparametric  elements  yielding  181  nodes  and 
147  elements.  The  properties  of  the  porous  medium  are  K  =  0.05  m/day 
and  effective  porosity  =  0.2.  Flow  vectors  for  steady  state  conditions 
are  depicted  in  Figure  1-9. 

The  length  of  arrows  represents  the  relative  magnitude  of  the  velocities 
It  is  interesting  to  note  that  most  of  the  water  is  drained  in  the  portion 
close  to  the  constant  head  zone.  Only  a  small  quantity  of  water  is  reach¬ 
ing  to  the  far  end  of  the  pipe. 
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INTERIM  MASS  FLUX  PROGRAM 


9.4.1  ABSTRACT: 

This  FORTRAN  program  will  plot  a  bar  or  line  graph  of  mass  fliix  on 
the  Tektronix  4051  screen  or  the  4‘662  plotter  for  user-provided  input 
parameters:  contaminant  information,  transect  definition,  geohydrologic 
information,  and  output  specifications. 

The  mass  flux  is  computed  from  the  formula 
0  =  KLtI  X  C  X  U 

where  0  =  mass  flux  across  the  transect  in  gms/day 

K  =  permeability  of  saturated  zone  in  feet/day 
L  =  saturated  length  in  feet 
t  =  saturated  flow  thickness  in  feet 
I  =  gradient  of  the  water  table 

C  =  contaminant  concentration  in  mg/liter  or  yg/liter 
U  =  units'  adjustment 

The  above  formula  is  essentially  Darcy's  law  for  the  flow  of  ground 
water. 

This  program  is  an  interim  version  of  the  mass  flux  calculation. 
Only  the  sampling  and  analysis  chemical  file  is  used  as  input  for 
programmatic  search.  All  geohydrologic  information  is  furnished  by 
the  user. 

9.4.2  PROCEDURE 

This  section  has  been  written  for  the  person  who  is  not  familiar 
with  a  digital  computer  but  is  motivated  to  use  the  computer  as  a  tool 
in  his  or  her  daily  job.  An  example  will  be  explained  step  by  step  to 
enable  the  user  to  gain  sufficient  knowledge  and  confidence  in  the  use 
of  this  program. 

Before  reading  further  it  is  necessary  that  the  user  understand 
how  to  log  on  the  Univac  1108  computer  as  described  in  sections  1.1 
through  1.3  of  the  IR  Data  Management  User's  Guide. 

Ensure  that  the  4631  hard  copy  unit  and  the  4662  plotter  are  turned 
on  and  the  four  binary  switches  on  the  plotter  are  set  to  "32A3." 

The  bold  type  below  will  indicate  the  computer's  response.  The 
prompting  symbol  ^  means  that  the  user  is  expected  to  furnish  input 
to  the  program  from  the  keyboard. 
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Example:  User  wishes  to  plot  the  mass  flux  of  contaminant 

DBCP  during  1978  across  a  transect  starting  at  Well 
23001  and  continuing  4000  feet  east.  User  determines 
that  there  should  be  eight  intervals  along  the  transect 
and  that  he  does  not  want  a  standard  concentration 
graph. 

The  user  should  study  Figure  4.1  on  the  next  page  in  order  to  see 
a  graphic  representation  of  the  mass  flux  formula. 
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Type  @ADD  IR*GEOLPRO.MF  and  press  [RETURN]  in  order  to  start 
the  program.  The  computer  will  respond  with: 


EHO  ERS. 

READY 

READY 

READY 

READY 

MAP  29R1  SL73R1  a5/'03x79  98:25:18 


The  computer  will  print  the  program  title  and  first  question: 


MASS  FLUX 

*  RUN"  PARAMETERS  t  19  JAN  79 

t  t  t  t  t  t  t  BASIC  INPUT 

ENTER  INSTALLATION  CODE  <2  LETTERS) 


Acceptable  responses:  RM  -  Rocky  Mountain  Arsenal 


Type  "RM"  and  press  [RETURN]. 
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The  computer  will  print  the  second  question: 


DO  YOU  WANT  TO 
POINT, 
OR  2 


SPECIFY  THE  TRANSECT 
DISTANCE,  ANGLE  <1> 
POINTS  <2>? 


ENTER  1  OR  2 


> 


There  are  two  methods  of  specifying  the  transect: 

•  point  1,  distance  from  1  at  bearing  angle  a  from  grid  north 


1 

•  point  1  and  point  2 


Either  method  allows  the  point  to  be  entered  as  site  type  and  site  ID  or 
in  state  planar  coordinates. 

For  the  example,  type  ”1"  and  press  [RETURN]  , 

The  computer  responds: 


> 


DO  YOU  WANT  TO  USE 


R  POINT  1,  in 

SITE  TYPE  +  SITE  ID 

nn  mnonTWOTPS 


<1) 


ENTER  1  OR  2 


Type  "1"  and  press  [RETURN]. 


( 
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The  computer  will  next  ask  for  site  type  and  site  ID. 


ENTER  SITE  TYPE 

examples: 


<4  LETTERS} 
BORE,  MELD 


> 


Acceptable  responses:  BORE 

WELL 

Type  "WELL"  and  press  [RETURN]. 


The  computer  responds : 


> 


ENTER  SITE  ID  <1  TO  10  LETTERS 
OR  DIGITS}  EXAMPLE  91001) 


Type  "23001"  and  press  [RETURN] . 


The  next  cwo  questions  that  the  computer  will  ask  concern  distance 
and  bearing  angle. 

ft 


ENTER  DISTANCE  FROM  POINT  1  IN  FEET 

<l  TO  5  DIGITS  +  DECIMAL  POINT} 
EXAMPLE;  5280.) 


Acceptable  response:  number  with  decimal  point 
Type  "4000."  and  press  [RETURN]. 

The  computer  responds: 


ENTER  BEARING  ANGLE  FROM  POINT  1  IN  DEGREES 

CLOCKWISE  FROM  NORTH}  0.  =  GRID  NORTH 

<1  TO  3  DIGITS  +  DECIMAL  POINT;  EXAMPLE:  99.) 


> 
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number  with  decimal  point. 


Acceptable  response: 

Type  "90."  and  press  [RETURN]. 


-EHIER  IHTERggL^HIDTH^«LOHG^THE^TR«HSECT 
EXAMPLE;  200.) 


IN  FEET 


> 


Acceptable  response:  number  with  decimal  point. 

This  number  is  W  in  Figure  4.1.  Type  "500."  and 
press  [RETURN]. 


Next  the  computer  will  ask: 


EKTER 


PERPENDICULAR  DISTANCE  FROM  THE  TRANSECT 
<1  TO  4  DIGITS  +  DECIMAL  POINTJ 

example:  200.) 


Acceptable  response:  number  with  decimal  point. 


IH  FEET 


This  number  is  P  in  Figure  4.1.  Type  ”1000.  and 
press  [RETURN].  Note  that  at  this  point  user  has 
built  *a  box  of  length  4000  feet  and  width  2000  feet 
comprised  of  8  smaller  boxes  as  shown  below  in 


Figure  4  .2 
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The  computer  will  compute  the  nxjmber  of  intervals  along  the 
transect: 


YOU  HAVE  8  INTERVALS. 

DO  YOU  WANT  TO  SPECIFY  A  NEW  INTERVAL  WIDTH?  <Y  OR  N) 


Type  "N"  and  press  [RETURN] . 

The  next  group  of  questions  concern  contaminant  input: 


>Y 


t  t  t  t  i  t  t  CONTAMINANT  INPUT 

DO  YOU  WANT  TO  SPECIFY  A  CHEMICAL 
TO  GET  A  COMPUTED  GRAPH? 


<Y  OR  H> 


Type  "Y"  and  press  [RETURN] . 


The  computer  responds: 


> 


ENTER  CHEMICAL  TEST  NAME  <1  TO  6  LETTERS 
OR  digits;  EXAMPLE  DBCP> 


Acceptable  response:  chemical  test  name 

This  contaminant  is  taken  from  section  2.1  of  the  IR  Data 
Management  User's  Guide.  Type  "DBCP"  and  press  [RETURN] . 


ENTER  SITE  TYPE  FOR  THE  SAMPLING  <4  LETTERS} 
EXAMPLES:  BORE,  WELL) 

> 


Acceptable  responses:  BORE 

WELL 

Type  "WELL"  and  press  [RETURN] . 
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c 


The  computer  will  ask  for  the  date  range: 


ENTER  BEGINNING  JULIAN  DATE  <0 

OR  5  digits;  example:  78244) 


Acceptable  response: 

where  YY  = 
DDD  = 
0  = 


YYDDD 

two-digit  year 
day  of  the  year 
earliest  date  on  file 


or  0 


Type  "78001"  and  press  [RETURN] . 
The  computer  responds: 


ENTER  ENDING  JULIAN  DATE  <0 
OR  5  DIGITS;  EXAMPLE: 


FOR  LATEST 
78244) 


c 


Acceptable  response:  YYDDD 


or  0 


■  where  YY  =  two-digit  year 
DDD  =  day  of  the  year 

0  =  latest  date  on  the  file 

Type  "78365"  and  press' ‘[RETURN]  . 

Next  the  computer  will  ask  how  to  deal  with  multiple  samples  for 
one  site  type  +  site  ID  that  arise  during  the  search  procedure: 


WHERE  MULTIPLE  SAMPLES  EXIST  AT  JHE  SUE  TYPE+ID 
AND  WITHIN  THE  DATE  RANGE  CHOSEN, 

DO  YOU  WANT  HIGHEST,  LATEST,  OR  MEAN?  <H,  L,  OR  M) 

NOTE'  USE  OF  MEAN  UALUES  OUER  LONG  INTERUALS 

OF  TIME  OR  DEPTH  HILL  MASK  TRENDS)  THEREFORE, 
PROCEED  WITH  CAUTION. 


) 


Acceptable  responses: 

H  -  Highest  sample  value  at  the  location  for  the  date 
range  specified 

L  -  Latest  (most  recent)  sample  value  at  that  location 
M  -  Mean  or  weighted  average  of  sample  values  at  that  location 
Type  "H"  and  press  [RETURN]. 
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At  this  point  the  user  can  specify  a  second  curve  to  be  plotted; 
i.e.,  the  standard  concentration  curve  on  the  same  set  of  x-  and  y-axes, 
which  can  be  useful  in  comparing  high  and  low  values  of  mass  flux.  To 
specify  standard  concentration,  user  will  have  to  enter  chemical  test 
name,  units,  and  value. 


> 


DO  YOU  WONT  TO- SPECIFY  A  CHEMICAL 

TO  GET  A  STANDARD  GRAPH?  <Y  OR  H> 


The  user  does  not  want  a  standard  graph,  so  type  "N"  and  press 
[RETURN] . 
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The  next  group  of  questions  concerns  geohydrologic  input; 


t  t  t  t  t  *  t  GEOHYDROLOGIC  INPUT 

*  *  *  *  t  *  *  FOR  IHTERUAL  1 

ENTER  PERMEABILITY  '  IN  FEET^DAY 

<1  TO  4  DIGITS  +  DECIMAL  POIHTJ 
EXAMPLE:  4.44) 


Acceptable  response:  nximber  with  decimal  point 
Type  "100.”  and  press  [RETURN]. 


ENTER  HATER  TABLE  GRADIENT 

(DECIMAL  POINT  +  1  TO  6  DIGITS? 


EXAMPLE: 


.080155) 


> 


Acceptable  response;  decimal  point  and  number 
Type  ".01"  and  press  [RETURN]. 


ENTER  SATURATED  THICKNESS  IN  FEET 

<1  TO  3  DIGITS  +  DECIMAL  POINT? 


example: 


90.  ) 


Acceptable  response:  number  with  a  decimal  point 
Type  "10."  and  press  [RETURN], 


> 


ENTER  SATURATED  LENGTH 

<1  TO  5  DIGITS  + 
example:  5280.) 


IN  FEET 

DECIMAL  point; 


Acceptable  response:  number  with  decimal  point 
This  number  does  not  have  to  equal  interval  width  U. 
Type  "500."  and  press  [RETURN]. 
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The  preceding  four  questions  will  be  repeated  for  each  interval  until 
the  user  has  entered  all  geohydrologic  information. 


t 


Next  the  computer  responds: 


‘t  Z  t  t  %  t  OUTPUT  SPECIFICATION* 

,0  vou  -gNT^VOUR^OUVUT^TO  BE  CRBBH 


Acceptable  responses: 

B  -  Bar  graph  for  each 
L  -  Line  graph  for  each 
Type  "L"  and  press  [RETURN]  . 


Next,  a  table  of  input  parameters  will  be  printed: 


*  RUN  PARAMETERS  t 


4  MAY  79 


IHS:  Rt1 

COHTAM 

INTERU  width: 


l:  OBCP  78001-78365 

500,  perpendicular: 


H 

1000. 


TRANSECT  LENGTH:  4000.  AT  90.  DEG 


INTERU 

1 

2 


PERM 

100.0 

100.0 


GRAD  THICK  LENGTH 
.0100  10.0  tttttt 

.0100  10.0  500. 


Assume  an  error  has  been  made  at  the  value  indicated  by  asterisks. 


DO  YOU  WANT  TO  CORRECT  ANY  GEOHYDROLOG IC  UALUES?  <Y  OR  H) 

> 


Type  "Y"  and  press  [RETURN] . 


> 


ENTER  IHTERUAL  NUMBER 

<2  DIGITS  INCLUDING 


LEADING  ZEROS  EXAMPLE: 


03) 


Type  "01"  and  press  return. 
Note  the  leading  zero. 
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ENTER  COLUMN  NUMBER 

<1  DIGIT;  1«PERM,2*GRAD,3=THICK,4=LENGTH) 

> 

Type  "4"  and  press  [RETURN] . 

EHTER  CORRECTED  ORLUE^  ^ 

> 

Now  the  user  should  enter  the  new  value;  type  ”500."  and  press  [RETURN]. 

DO  YOU  WANT  TO  CORRECT  ANY  GEOHYDROLOGIC  VALUES?  <Y  OR  H> 


Type  "N"  and  press  [RETURN] .  The  new- table  of  input  parameters 
will  be  printed. 


t  RUN  PARPMETER?  t  8  MAY  79 
IHS*  Rn 

COHTAM  1;  DBCP  78001-78365  H 
IHTERU  WIDTH;  580.  PERPENDICULAR;  1000. 
TRANSECT  LENGTH;  4000.  AT  90.  DEG 


IN7ERV 

PERM 

GRAD 

THICK 

LENGTH 

1 

108.0 

.0100 

10.0 

500. 

2 

100.0 

.0100 

10.0 

500. 

3 

100.8 

.•0100 

10.0 

500. 

4 

100.8 

.0108 

10.0 

508. 

5 

100.0  ■ 

•.0100 

10.0 

500. 

6 

100.0 

.0100 

18.0 

500. 

7 

100.0 

-.0100 

10.0 

500. 

8 

100.0 

.8100 

10.0 

580. 

100 


The  computer  will  search  the  IR  Data  Base  for  candidate  subschema 
records. 

The  resultant  graph  will  be  drawn  on  the  screen  and  the  program  will 
stop.  If  the  user  desires  a  hard  copy,  press  [MAKE  COPY];  then  press 
[RETURN]  to  continue. 


The  computer  will  then  ask  if  the  same  graph  should  be  plotted  on 
the  4662  plotter: 

DO  YOU  HftHT  THIS  GRAPH  OH  THE  4662  PLOTTER?  <Y  OR  N> 

> 

Type  "Y"  and  press  [RETURN] . 

LOAD  1  SHEET  OF  PAPER}  THEN  PRESS  CRETURH3 

>  “  ■ 

The  computer  will  pause  so  that  the  user  may  load  a  sheet  of  plotting 
paper.  To  continue,  press  [RETURN]  and  then  the  graph  will  be  plotted. 

After  the  plot  the  computer  will  ask  if  the  user  wishes  to  stop. 

If  not,  user  should  type  "N"  and  answer  the  same  questions  as  before. 

When  user  has  completed  all  his  computer  runs,  type  (3ADD  IR*GEOLPRO.EXIT 
and  press  [RETUTIN]  to  clear,  the  computer;  the  computer  will  respond: 


EXITS 

FURPUR  27R3  E33  SL73R1  05/04^79  11:18:35 
EHD  ERS. 
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Note: 


If  there  are  problems  with  the  program, 
by  typing  §§X  TIO  and  pressing  RETURN. 


stop  execution 
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COMPUTED  MASS  FLUX  OF  DBCP  AT  ROCKY  MT  ARS 
FROM  78001  TO  78365 


9.4.3 


SAMPLE  OUTPUT 


Note:  Computed  mass  flux  -  solid  line 


Standard  mass  flux 


dashed  line 


CD 


U. 


00  KS) 
00 

wo 

U  il 
>:> 


^  szcLfftfji  U.-J3X  L3  2:  \  ^:::4<x>'»<^ouooo 
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COMPUTED  MASS  FLUX  OF  DBCP  AT  ROCKY  MT  ARS 
FROM  78001  TO  78363 


C 


I- 

u. 


C\J^ 
03  VO 
rt  rw 
-*-1  ID 

coc^ 

CUCD 
U  R 
30 


2:<xco  w  u. 

L, 


=330 


C32:N^<3:3-^*C3li30Q3 
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COMPUTED  MASS  FLUX  OF  DBCP  AT  ROCKY  MT  ARS 
FROM  78001  TO  78365 

STANDARD  MASS  FLUX  OF  DBCP  AT  ROCKY  MT  ARS 
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X»2181482 

Y=0 195761  LENGTH  ALONG  TRANSECT  <FT) 


Can|)bell/inas/2041 
14  Nov  79 


KQCIH-IS 

SUBJECT:  Proposed  Position  on  Fluoride  Txeataent  at  RMA 


Consnaider 

Rodey  Mountain  i\rsenal 
Canasice  City,  00  80022 


1.  Reference  meeting  at  RM^,  4  Oct  79,  between  representatives  of 
FMA,  WES,  and  this  .^ency,  subject  as  above. 

2.  At  re&renced  meeting  results  of  recently  coo^leted  gedfaydrologie 
survey  tasks  at  the  north  boiaidaiy  of  IMA.  and  results  of  fluori<te 
rsnioval  pilot  tests  were  discussed.  Because  of  the  hi^  costs  of 
fluoride  treatment  and  the  uncertainty  associated  with  the  requirenient 
to  treat  for  fluoride  t^xm  systm  expansion,  a  decision  was  reached  / 
to  advise  the  State  of  Colorado  that  the  US  ktesf  does  not  plan  to 
construct  a  fttoride  removal  systaa  at  this  time.  The  ongoing  design 
of  a  fluoride  treatment  process  will  be  continued  to  allow  rapid 
isplem^tatioa  if  needed  at  a  latter  date. 

3.  Per  agre^nsnts  smde  with  !9'!A,  attadied  at  inclosure  1  is  a 
suggested  letter  for  transsdttal  to  Colora^  Departmmt  of  Health. 
Inclosure  2  is  sipportive  data  compiled  by  USA3HAMA  on  the  expected 
fluoride  contaminant  loading  to  the  ejq>anded  north  boundary  control 
systen. 

4.  Reciuost  this  informaticsi  be  provided  to  the  State  of  Colorado 
as  soon  as  possible. 


2  Incl 
as 


FRANK  A.  JONES,  JR. 
Colonel  ,  QnlC 
Cesananding 


DATE  DUE  j 

